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Section  1 


(INTRODUCTION  AND  SUMMARY 

The  accurate  prediction  of  the  fallout  hazard  that  results 
from  a subsurface,  surface,  or  low  air  burst  is  perhaps  one  of  the 
most  complex  tasks  that  has  faced  military  analysts.  It  encompasses 
many  fields  of  technology  ranging  from  radio -chemistry  and  thermo- 
dynamics to  meteorology  both  micro-scale  and  synoptic.  However, 

5 

the  mere  fact  that  a component  of  the  problem  is  difficult  to  under- 
stand and  model  does  not  necessarily  mean  that  progression  from  a 
simple  algorithm  to  a sophisticated  computation  changes  the  ultimate 
answer  significantly.  Hence,  there  is  a continuing  need  for  sensitivity 
analyses  to  indicate  both  how  accurate  a present  code  is  and  perhaps 
what  might  be  important  areas  in  need  of  improvement. 

This  report  presents  the  results  of  a study  on  the  sensitivity 
of  fallout  predictions  to  uncertainties  associated  first  with  the  burst 
parameters  (Section  3)  and  secondly  with  assumptions  and  models 
contained  in  the  overall  prediction.  The  basis  for  the  study  was  the 
calculation  of  the  DELFIC  code  (Department  of  Defense  Land  Fallout 
Prediction  System).  As  will  be  discussed  in  Section  2,  the  Ballistics 
Research  Laboratory  made  a series  of  runs  with  DELFIC.  Science 
Applications,  Inc.,  (SAI)  with  partial  modifications,  exercised  the 
Particle  Activity  and  Output  Processor  Modutes. 

In  addition  to  uncertainties  in  burst  parameters,  several 
other  areas  were  investigated.  The  first  of  these  is  the  size  distri- 
bution (Section  4)  of  the  entrained  dust  for  which  different  types  of 
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distributions  were  analyzed.  The  effect  of  variations  in  heights  of 
burst  is  also  discussed.  Section  5 deals  with  several  aspects  of  the 
rise  phase  including  time -dependent  cloud  dimensions  and  the  internal 
distribution  of  dust.  Section  6 is  devoted  to  the  Particle  Activity  Module 
with  emphasis  on  fractionation.  Finally,  Section  7 reports  on  analysis  of 
neutron-induced  soil  activation  which  may  be  significant  under  certain 
circumstances.  The  various  appendices  deal  with  several  problems 
and  inconsistencies  found  in  DELFIC  and  provide  a listing  of  a 
modified  output  processor. 


1. 1 SUMMARY 

This  section  presents  some  important  conclusions  of  the 

study. 


a.  Variations  in  the  analytic  form  of  the  size  distri- 
bution can  result  in  changes  in  close-in  dose  of 

a factor  of  20  and  a factor  of  5 for  long-range 
dose  levels  (Section  4. 1). 

b.  Present  HOB  algorithms  overestimate  dose 
levels  close-in  particularly  for  low  yields 
(Section  4. 2). 

c.  DELFIC  predicts  a larger  and  less  dense  cloud 
than  does  the  analytic  flow  field  model,  VORDUM, 
and  does  not  loft  particles  greater  than  1 cm  (at 
300  kt).  This  plus  the  differences  in  the  spatial 
distribution  of  the  dust  between  VORDUM  and 
DELFIC  result  in  a discrepancy  in  the  predicted 
dose  of  approximately  a factor  of  2 and  in  the  pre- 
dicted time  cf  arrival  of  about  15  minutes  (Section  5.1). 

d.  The  differences  in  rise  parameters  between  DELFIC 
and  VORDUM  (with  semi-empirical  equations)  result 
in  slightly  under  an  order -of -magnitude  discrepancy 
in  the  predicted  dose  rate  for  the  30  kt  case.  At 
higher  yields,  there  is  better  agreement  (Section  5.2). 
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e.  An  alternative  formalism  of  the  Freiling 
fractionation  model  changes  the  dose  by 
about  50  percent  (Section  6. 1). 

f.  Assuming  the  debris  to  be  purely  refractory 
for  all  mass  chains  satisfactorily  reproduces 
the  more  complex  Freiling  model,  at  least  for 
a power  law  size  distribution  (Section  6. 2). 

g.  The  neutron  energy  spectra  is  an  important 
parameter  in  determining  the  fraction  of  neutrons 
captured  in  crater  material  that  is  later  entrained 
into  the  clocd  (Section  7. 1).  The  neutrons  above 
4 Mev  are  twice  as  efficient  in  contributing  to  the 
total  activity  of  entrained  dust  as  those  below 

4 Mev  (Section  7. 1). 

h.  A factor  of  5 variation  was  noted  in  the  induced 
activity  between  the  various  soil  types  examined 
(Section  7.2), 

i.  DELFIC  overestimates  the  contribution  of 
induced  activity  for  heights  of  burst  greater 
than  a few  feet  (scaled  to  1 kt)  (Section  7.4). 

j.  DELFIC  assumes  the  activated  soil  to  be 
refractory  in  the  fractionation  scheme.  If  one 
assumes,  with  some  justification,  that  it  is 
volatile  then  for  the  burst  configuration  examined, 
the  result  is  an  order-of -magnitude  decrease  in 
the  very  close-in  dose  levels  due  to  the  activated 
soil  component  of  the  fallout  (Section  7.0). 


f 
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Section  2 
PROCEDURE 

The  basic  DELFIC  runs  were  performed  at  the  Ballistic 
Research  Laboratories  (BRL)  and  material  in  the  form  of  data  tapes 
and  output  were  sent  to  SAI  for  further  computation  and  analysis.  A 
total  of  ten  yields  was  included.  These  were  30  kt,  100  kt, 

270  kt,  300  kt,  330  kt,  1 Mt,  2. 7 Mt,  3 Mt,  3. 3 Mt,  and  10  Mt.  The 
topographical  and  wind-field  data  were  identical  for  all  runs;  the 
only  input  differences  were  yield-dependent  parameters.  Additionally, 
the  300  kt  run  was  executed  using  the  horizontal  subdivision  option 
in  the  cloud  rise  module;  all  other  runs  were  performed  with  no 
horizontal  cloud  subdividing.  Table  2. 1 contains  the  invariant  input 
parameters  for  all  runs,  while  Table  2. 2 contains  the  variant 
parameters. 

A complete  input  deck  for  DELFIC  contains  many  more 
parameters,  some  involving  machine-related  variables  such  as 
number  of  tape  drives  available,  and  others  being  standard  variables 
such  as  atmospheric  profiles.  Since  DELFIC  is  a growing  program 
in  the  sense  that  its  various  modules  are  being  continuously  updated, 
no  readily  available  operator’s  manual  exists.  The  personnel  at 
BRL,  on  consultation  with  SAI,  provided  the  necessary  input  variables 
not  listed  in  the  following  tables. 

A critical  input  parameter  is  the  specification  of  particle 
size  distribution.  Traditionally,  users  of  DELFIC  have  assumed 
a log-normal  distribution.  However,  as  will  be  discussed  in  Section  4, 
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Table  2. 1.  Invariant  DELFIC  Input  Parameters 


Number  of  Particle  Size  Classes 

40 

Type  of  Size  Distribution 

a"3,  (power  law) 

Soil  Specific  Gravity 

2.5 

Soil  Solidification  Temperature 

1673  °K 

Height  of  Burst 

0.0 

Type  of  Atmosphere 

Standard  U.  S 

. -45°  N July 

Time  Range  of  Interest 

0 - 48  Hr 

Type  of  Topography 

Planar  at  MSL 

Fission  Type 

U235  HE 

Capture  to  Fission  Ratio 

0.0 

Ground  Roughness  Factor 

.5 

Wind  Hodograph  Extending  Over 

Complete  Range  of  Interest 

Winds  for  60N  - 20E  Winter 

Vector  Altitude,  ZV(J)  (m) 

VX(J)  (m/sec) 

VY(J)  (m/sec) 

0.00000 

4. 55700+00 

-8. 19000-01 

3.04800+03 

4.55700+00 

-8. 19000-01 

6.09600+03 

6.75800+00 

-2.45800+00 

9. 14400+03 

1.00860+01 

-1.79200+00 

1.21920+04 

1.20320+01 

-4.40300+00 

1.52400+04 

1.20320+01 

-4.40300+00 

1.82880+04 

1. 51040+01 

-2. 66200+00 

2.13360+04 

1.74080+01 

-3. 12300+00 

2.43840+04 

2.04800+01 

0.00000 

2.74320+04 

2. 30400+01 

0.00000 

3.04800+04 

2.56000+01 

0.00000 
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Table  2. 2.  Variant  DELFIC  Input  Parameters 


Yield  (kt) 

CAY 

FW  (kt) 

KDI 

EMITN 

IRAD 

30 

7.28xl07 

30 

m 

mm 

0 

100 

2.11x10® 

100 

is. 

0 

270 

5.29x10® 

270 

40 

1.5 

0 

300 

5.29x10® 

300 

40 

3 

330 

5. 29x10® 

330 

40 

1 1 

o 

1000 

1.61x10® 

500 

40 

H 

mm 

2700 

4.25x10® 

1350 

/ 

n 

3000 

4.25x10® 

1500 

B] 

m 

3300 

4.25x10® 

1650 

X9 

Wm 

10000 

1.23x10*® 

5000 

50 

■n 

CAY  --  Mass  Normalization  Factor 

FW  = Fission  Yield 

KDI  = Number  of  Wafers  in  Cloud  Rise  Module  Per  Size 
Class 

EMITN  = Number  of  Neutrons  Per  Fission 
IRAD  = Horizontal  Wafer  Subdivision  Factor 
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there  is  conclusive  evidence  that  the  distribution  !ar  particles  with 
diameters  above  a lew  hundred  microns  obeys  a power  law.  Con- 
sequently, the  base-line  calculations  generated  by  BRL,  were  per- 
formed with  the  distribution  on  particle  size  of  the  following  form: 


dn 

la 


where  a is  diameter  and  dn/da  is  the  number  of  particles  per  unit 
diameter.  The  normalization  to  the  total  mass  burden,  M,  of  the 
cloud  is  made  through  the  input  variable  CAY  as  follows: 


where  p is  the  bulk  density  of  the  soil  and  the  particles  are  assumed 
to  be  spherical.  The  limits  on  the  integration  were  taken  to  be  0 
and  10 V-  The  total  mass  is  expressed  in  grams.  Table  2. 2 lists 
the  values  of  CAY  for  the  various  yields. 

The  first  case  exercised  was  for  a yield  of  300  kt.  It  was 
found  that  the  DELFIC  formalism  for  gravitational  settling  was 
invalid  for  particles  greater  than  1 cm.  For  yields  in  the  kiloton 
and  megaton  range,  particles  of  this  size  and  larger  are  expected 
to  be  entrained.  Using  the  VORDUM  code^\  which  describes  the 
velocity  field  in  and  around  the  rising  fireball,  Figure  2. 1 was  gen- 
erated for  the  largest  particle  entrained  versus  yield  (for  surface 

bursts  only).  DELFIC  (in  the  RSXP  subroutine)  makes  use  of  an 

(2) 

approximate  expression  developed  by  Davies'  for  the  drag  coefficient 
as  a function  of  Reynold's  number  that  does  not  apply  to  large  par- 
ticles (greater  than  1 cm).  Because  of  this,  the  value  of  CAY  for  all 
subsequent  runs  was  chosen  so  that  the  upper  limit  on  particle  size  was  1 cm. 
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Figure  2. 1.  Largest  Particle  in  Main  Cloud  as  Predicted  with  VORDUM 


The  data  sent  to  SAI  for  analysis  and  further  computation 
were  in  the  form  of  magnetic  tapes  and  computer  printouts.  The 
data  on  the  tapes  were  the  stabilized  cloud  description  from  the 
Cloud  Rise  Module  (option  KRX  = 1),  and  the  grounded  particles  tape 
(tape  IPOUT).  Initially,  tape  difficulties  due  to  machine  incon- 
sistencies arose,  but  all  the  necessary  data  were  finally  transmitted. 

Additionally,  a tape  containing  card  images  of  the  complete 
DELFIC  code  was  brought  to  SAI.  BRL  uses  a UNIVAC  1108,  while 
SAI  has  access  to  a CDC  6400,  so  the  DELFIC  code  had  to  be 
modified  to  take  into  account  the  different  FORTRAN  compilers.  The 
Output  Processor  Module,  with  the  associated  PAM  routines,  was 
successfully  run  at  SAI  using  the  appropriate  grounded  particle 
tapes.  The  output  produced  was  identical  to  that  supplied  by  BRL, 
so  the  transfer  of  the  Output  Processor  was  successful. 

The  computer  printout  was  also  used  to  identify  an  appro- 
priate ''hotline"  for  the  wind  field  used.  It  appeared  that  this 
"hotline"  occurred  along  a line  345°  to  the  positive  x-axis  for  the 
larger  yields.  A line  of  350°  was  used  for  the  30  kt  case,  and  a 
line  of  347°  was  used  for  the  100  kt  case. 

The  form  of  the  standard  DELFIC  output,  a rectangular 
map,  was  inappropriate  for  purposes  of  a sensitivity  study  along  the 
hotline  and  along  lines  parallel  to  the  hotline  but  removed  from 
ground  zero  by  a stabilized  cloud  radius.  For  this  reason,  a program 
was  written  to  produce  the  desired  output  along  any  specified  line. 

This  program  '.s  called  LINTAP,  and  a listing  of  it  appears  in 
Appendix  C.  LINTAP  uses  all  the  associated  PAM  subroutines,  as 
well  as  subroutine  GETSET  from  the  Output  Processor.  The 
functions  that  were  performed  by  the  rest  of  the  Output  Processor 
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are  performed  in  subroutine  LINOUT.  The  inputs  to  LINTAP  are 
the  grounded  particles  tape  (IPOUT),  as  well  as  the  inputs  to  the 
PAM  module  and  a modified  version  of  the  inputs  to  the  standard 
Output  Processor.  Table  2. 3 is  an  explanation  of  the  inputs  to 
LINTAP,  while  Table  2. 4 is  an  explanation  of  the  various  requests 
(variable  NREQ).  It  was  desired  to  keep  modifications  of  the  Output 
Processor  to  a minimum,  which  is  the  reason  for  the  similarity  of 
the  input  deck  structures.  The  output  of  LINTAP  is  a tape  suitable 
for  input  to  a CALCOMP  program,  if  requested  (LDIS  4 0),  and/or 
the  values  of  any  of  the  16  quantities  along  the  specified  line.  If  a 
tape  is  requested,  the  information  is  written  on  Tape  4. 

The  logic  of  LINTAP  is  as  follows.  The  grounded  particles 
tape  contains  the  X and  Y coordinates  of  the  center  of  the  wafer  when 
it  hits  the  ground,  as  well  as  the  particle  diameter  associated  with 
the  wafer,  the  time  of  impact,  the  mass  per  unit  area,  and  a number 
(IRBZ)  from  which  the  radius  of  the  wafer  can  be  calculated.  The 
standard  Output  Processor  assumes  rectangular  wafers,  where  the 
Cloud  Rise  Module  assumes  cylindrical  wafers.  This  is  handled 
within  the  Output  Processor  by  adjusting  the  wafer  radii  to  conserve 
area  (hence  mass).  In  the  standard  Output  Processor,  a value  of  886 
for  the  variable  IRBZ  corresponds  to  a wafer  which  is  as  wide  as  the 
visible  cloud.  This  value  of  886  c—responds  to  (jt/4)2  x 1000,  the 
diameter  correction  factor  between  circles  and  squares  containing 
the  same  area.  LINTAP  assumes  the  wafers  are  cylinders,  and  the 
associated  IRBZs  are  adjusted  accordingly. 

Figure  2. 2 illustrates  the  actual  program  logic.  The 
specification  of  XMIN,  YMIN,  ALPHA,  RANMAX,  and  DGX  totally 
specifies  the  line  and  points  1 through  24  on  the  line.  The  program 
then  checks  each  wafer  for  intersection  points.  Wafer  #1  does  not 
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Table  2.3.  Inputs  to  LINTAP 
Tapes:  Tape  3 = IPOUT 

Card  Description  Format 

1 Run  Identifier  12A6 

2 Blank 

3*  IC(I),  1=1,  18  1814 

IC(17)  > 0 Stops  after  printing  IPOUT 
IC(17)  = 0 Proceed  with  Job 
IC(18)  > 0 Print  IPOUT 
IC(18)  = 0 Do  not  print  IPOUT 


DIFCON  F10.3 


5 IH,  IV  1814 


Run  Identifier  12A6 

Blank 


8 Blank 

9*  CAPFIS,  EMITN  2F10.3 

10*  FISSID  A6 

11  LLL,  LDIS  214 


LLL  = No.  of  lines  desired 

LDIS  > 0 Write  tape,  print  output 

LDIS  = 0 Do  not  write  tape, 
print  output 

LDIS  < 0 Write  tape,  do  not  print 
output 

(Continued) 
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Table  2. 3.  Inputs  to  LINTAP  (Continued) 


Tapes:  Tape  3 = IPOUT 

Card 

Description 

Format 

12* 

XMIN,  YMIN  ALPHA,  RANMAX,  DGX, 
GRUFF 

6F10.4 

XMIN,  YMIN  = Starting  point  of 
line  (m) 

ALPHA 

= Angle  of  line  to 
horizontal  (°) 

RANMAX 

= Maximum  range 
desired  (m) 

DGX 

= Increment  along 
line  (m) 

GRUFF* 

= Ground  roughness 
factor 

13* 

KKK1 

14 

•v 

+ 

KKK1 

= Number  of  output 
requests  per  line 

NREQ,  Tl,  T2, 

MASCHN 

14, 

2F10.3, 

14 

NOTE:  Cards  1 through  10  are  identical  to  the  cards  required  for 
the  standard  Output  Processor. 

* Signifies  input  in  the  first  ten  cards  that  is  still 
utilized  in  LINTAP. 

Signifies  cards  within  a loop  defined  by  LLL,  i.  e. , 
LINTAP  requires  LLL  sets  of  these  cards. 

+ Signifies  cards  within  a loop  defined  by  KKK1,  i.  e. , 
LINTAP  requires  KKK1  cards  14. 
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Table  2.4.  Available  Computation  Codes 


Computation  Code 
NREQ 

Computation  Type 
Description 

1 

Count  of  wafers  covering  each  output 
point 

2 

Exposure  rate  normalized  to  time  H + 1 
hour 

3 

Exposure  rate  at  time  H + T1  hours 

4 

Integrated  exposure,  H + T1  to  infinity 
accounting  for  time  of  arrival 

5 

Integrated  exposure,  H + T1  to  H + T2 
accounting  for  time  of  arrival 

b 

Total  mass  deposited 

7 

Total  mass  deposited  from  time  H + T1 
to  H + T2 

8 

Integrated  exposure,  H + T1  to  H + T2 
assuming  all  particles  have  arrived  by 
H + T1  hours 

9 

Same  as  8 integrated  to  infinity 

10 

Concentration  of  an  individual  mass  chain 
(curies/m^) 

11 

Time  of  arrival  of  first  fallout  particle 

12 

Time  of  deposit  of  last  fallout  particle 

13 

Smallest  particle  size  deposited 

14 

Largest  particle  size  deposited 

15 

Mass  from  particles  in  size  range  T1 
to  T2 

16 

H + 1 hour  "normalized”  exposure  rate 
resulting  from  particles  in  size  range 
T1  to  T2  microns 
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intersect  and  is  ignored  while  wafers  #2  and  #3  do  intersect. 
Following  the  output  processor  logic,  points  9 through  13  are  incre- 
mented equally  according  to  the  mass  (or  time,  etc. ) associated 
with  wafer  #2,  and  points  12  through  18  are  incremented  equally 
according  to  the  mass  associated  with  wafer  #3.  This  process  is 
continued  until  all  wafers  have  been  tested,  and  the  results,  values 
at  points  1 through  24,  are  then  printed.  If  another  request  exists, 
the  program  rewinds  the  grounded  particles  tape  and  starts  again. 

The  results  from  program  LINTAP  were  compared  with  the 
maps  produced  by  the  Output  Processor,  and  although  a direct  one- 
to-one  comparison  is  not  possible  (because  of  the  above-mentioned 
squares -circles  logic),  the  differences  are  insignificant.  LINTAP 
was  used  in  all  the  sensitivity  analyses. 

When  integrated  dose  (NREQ  = 4 or  5)  is  requested  of  the 
Output  Processor,  the  run  almost  always  requires  astronomical 
amounts  of  time.  For  example,  a grounded  particles  tape  of  34 
particles  needs  118  sec  of  CDC  6400  executive  time.  Most  of  the 
tapes  used  in  this  study  contained  approximately  2000  particles,  so 
a faster  method  of  calculation  was  developed  for  LINTAP.  This 
long  time  is  needed  for  integrated  dose  rates  because  each  particle 
■ requires  a call  to  PAM2  and  its  associated  routines,  while  for  all 

| other  requests,  PAM2  is  called  once,  if  at  all,  and  the  data  stored 

? in  array  FP  are  used  for  all  particles.  In  effect,  the  integrated 

doses  are  the  only  time-of -impact  dependent  requests.  A modifica- 
tion to  the  Output  Processor  logic  was  made  within  LINTAP  to 
increase  the  speed  of  execution  of  integrated  dose  requests,  with 
little  loss  of  accuracy. 
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PAM2  calculates  the  integrated  dose  for  all  particle  sizes 
from  time  T1  to  time  T2  (either  infinity  or  a specific  value  depending 
on  NREQ  = 4 or  5),  and  stores  this  information  in  array  FP.  The 
Output  Processor  s^ts  T1  equal  to  the  time  of  impact  of  the  specific 
particle,  calculates  the  integrated  doses  for  all  particle  sizes,  and 
then  just  uses  the  one  value  for  the  specific  particle  size.  LINTAP 
calculates  the  integrated  doses  for  all  particle  sizes  at  a number  of 
specific  times,  and  saves  these  values  in  array  FCAN.  The  time  of 
impact  is  then  used  as  the  parameter  to  interpolate  between  the 
various  times  without  having  to  resort  to  PAM2.  Figure  2.3  is  a 
simplified  sketch  showing  this  process. 

The  times,  in  hours,  at  which  the  dose  is  calculated  are  0, 

0.25,  0.5,  0.75,  1,  2,  4,  6,  etc.,  in  two-hour  increments  until  time 

of  cessation  or  time  requested  (NREQ  = 5)  is  reached.  These  times 

were  chosen  on  the  basis  of  the  experimental  evidence  which  showed 

-1  2 

that  the  dose  from  early  fallout  followed  a T ’ dependence,  and 
from  actual  LINTAP  development  runs  which  showed  that  better  than 
two-hour  resolutions  were  needed  for  early  times.  Also,  on  the  basis 
of  the  above-mentioned  power-law  dependence,  the  interpolation  is 
done  as  a log -log  procedure. 

This  process  was  checked  using  the  3 Mt  grounded  particles 
tape.  Using  the  output  received  from  BRL,  regions  of  the  map  were 
identified  which  contained  very  few  particles.  A version  of  LINTAP 
was  used  which  calculated  the  integrated  dose  from  time  of  arrival 
to  48  hours  after  burst  (NREQ  = 5)  according  to  the  Output  Processor 
procedure.  The  dose  from  these  regions  was  also  calculated  using 
the  above  procedure.  Table  2. 5 summarizes  the  results  of  these 
comparisons.  It  can  be  seen  that  the  doses  are  predicted  by  the  Output 


28 


Three  Particle  Size  Classes 
Particle  A of  diameter  #1  impacts  at  T = 1.5  Hr. 
Particle  B of  diameter  #2  impacts  at  T = 0. 7 Hr. 
Particle  C of  diameter  #3  impacts  at  T = 3.0  Hr. 
Particle  D of  diameter  #2  impacts  at  T = 2. 5 Hr. 

Etc. 


♦Note:  Calculations  are  log-log  within  LINTAP, 
semi -log  is  just  used  for  demonstration 
purposes. 


5 6 

Time  (Hr) 


Figure  2. 3.  Simplified  Dose  Calculation  Within  LINTAP 
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Processor  over  a wide  range  of  impact  times.  Further  comparisons 
were  not  made  because  of  the  cost  involved  with  generating  the 
"standard"  doses. 


Table  2. 5.  Integrated  Dose  Calculation  Comparison 


Number  of 

Time  of  Impact 

Standard  Dose 

LINTAP  Dose 

Particles 

(Sec) 

(R) 

(R) 

l 

66970 

.28625186 

.28552985 

2 

2382  - 3689 

5. 53905285 

5.524743 
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Section  3 

SENSITIVITY  TO  BURST  PARAMETERS 

3.1  OBJECTIVE 

Even  in  the  most  elaborately  fielded  nuclear  test,  there  is 
some  uncertainty  in  the  total  yield  estimate,  as  well  as  the  total 
fission  and  fusion  yields.  In  an  operational  situation,  these  uncer- 
tainties in  total  yield,  fission  yield,  and  also  fission  type  could  be 
greatly  magnified.  This  section  investigates  the  effects  of  these 
uncertainties  on  the  predicted  fallout. 

3. 2 BASELINE  CALCULATIONS 

r The  dose  rate  and  total  dose  for  five  yields  spanning  the 

range  30  kt  to  10  Mt  were  determined  along  the  downwind  radial 
from  ground  zero,  (the  "hotline"),  and  along  two  lines  parallel  to  the 

i 

} hotline  but  displaced  from  it  by  the  visible  cloud  radius  at  stabilization. 

The  five  yields  were  30  kt,  100  kt,  300  kt,  3 Mt,  and  10  Mt. 

j Figures  3. 1 and  3. 2 are  the  dose  rate  normalized  to  H+ 1 hr 

* and  the  total  integrated  dose  48  hours  after  burst,  respectively,  for 

the  30  kt  device  along  the  three  lines.  In  these  and  subsequent  figures 
in  this  subsection,  the  solid  line  corresponds  to  the  hotline  output,  the 
'Vs"  correspond  to  the  parallel  line  displaced  to  the  left  of  the  hotline 
by  the  cloud  radius  at  stabilization,  and  the  "x's"  correspond  to  the 
line  one  stabilized  cloud  radius  to  the  right  of  the  hotline.  The  run 
parameters  listed  in  the  preceding  section  were  used  for  all  the 
yields  considered.  (Note:  The  irregularity  associated  with  all  the 
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LINTAP  output  is  due  to  the  division  of  the  cloud  into  a finite  number 
of  wafers. ) 

Figures  3.3  and  3.4  are  the  dose  rate  and  total  integrated 
dose  48  hours  after  burst,  respectively,  for  the  100  kt  device  along 
the  three  lines.  Similarly,  Figures  3. 5 and  3. 6 present  the  same 
information  for  300  kt,  while  Figures  3.7  and  3.8  correspond  to 
3 Mt  and  Figures  3. 9 and  3. 10  correspond  to  10  Mt. 

As  expected,  the  above  figures  show  that  the  downwind  pro- 
jection of  the  cloud  outline,  the  cloud  "shadow, " is  the  only  area  that 
need  be  considered.  Figure  3. 11  is  a plot  of  the  maximum  downwind 
extent  of  two  dose  levels  48  hours  after  burst  as  a function  of  yield. 
The  dose  levels  considered  are  150  R (the  mid-burdening  dose),  and 
450R  (the  mid-lethal  dose).  For  these  two  levels,  at  least,  the 
maximum  downwind  extent  is  not  a linear  function  of  yield,  which 
means,  for  example,  that  doubling  the  total  yield  will  not  double  the 
downwind  extent  of  these  levels.  Rather  the  range  for  either  dose 

p 

level  is  proportional  to  W where  W is  the  yield  and  P is  approxi- 
mately 0. 22. 

3. 3 TOTAL  YIELD  AND  FISSION  FRACTION  SENSITIVITY 

This  subsection  investigates  the  fallout  variations  due  to 
uncertainties  in  the  total  yield  and  fission  yield  on  the  order  of 
*10  percent  which  are  the  approximate  variations  that  could  be  ex- 
pected when  the  designed  yield  of  the  device  is  known.  The  two 
specific  yields  investigated  were  300  kt  and  3 Mt. 

The  first  effect  studied  was  total  yield  sensitivity. 

Figure  3. 12  is  a plot  c'  ose  rate  along  the  hotline  normalized  to 
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Figure  3.3.  Normalized  Dose  Rate  Along  Angle  (Deg)  « 347;  Yield  ■ 100  kt 
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Figure  3. 6.  Integrated  Dose  Along  Angle 
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Figure  3.7.  Normalized  Dose  Rate  Along  Angle  (Deg)  = 345;  Yield  = 3 Mt 


tgure  3.8, 


Figure  3. 9.  Normalized  Dose  Rate  Along  Angle  (Deg)  * 345;  Yield  * 10  Mt 


Figure  3. 10.  Integrated  Dose  Along  Angle  (Deg)  = 345;  Yield  * 10  Mt 


Variations  in  Predicted  Dose  Rate  Due  to  *10  Percent 
Uncertainties  in  Total  Yield 


H+l  hour  for  three  yields,  2.7  Mt,  3.0  Mt,  and  3.3  Mt,  where  both 
the  fission  yield  and  total  yield  have  been  varied  proportionally 
(*10  percent  around  the  3 Mt  baseline  case).  Figure  3.13  Is  the 
hotline  total  integrated  dose  48  hours  after  burst  for  the  same  three 
yields.  The  doses  and  dose  rates  for  the  three  devices  are  tabeled 
by  yield  on  the  plot  itself.  It  can  be  seen  from  these  figures  that  the 
close-in  fallout  at  any  specific  range  varies  at  most  by  *10  percent 
with  yield,  and  that  the  long-range  fallout  appears  to  vary  by  a 
greater  amount,  but  that  this  latter  effect  is  most  probably  caused 
by  the  limited  number  of  cloud  wafers  used  in  the  DELFIC  formulism, 
and  thus  not  be  greater  than  the  close-in  variation  if  a greater  number 
of  small  particle  wafers  were  used. 

Similarly,  Figures  3. 14  and  3.  IS  are  the  respective  H+l  hr 
normalized  hotline  dose  rates  and  total  integrated  doses  48  hours 
after  burst  for  the  270  kt,  300  kt,  and  330  kt  devices.  The  variations 
of  the  close-in  fallout  are  again  within  *10  percent,  and  the  same 
conclusions  regarding  the  long-range  fallout  can  also  be  made. 

An  important  conclusion  regarding  the  future  execution  of 
the  DELFIC  code  can  be  made  from  these  figures.  To  accurately 
describe  the  long-range  fallout,  many  more  small  particle  wafers 
should  be  included  in  the  run.  Perhaps  an  option  of  specifying  the 
number  of  wafers  as  a function  of  particle  size  could  be  incorporated 
in  the  code. 

The  next  effect  investigated  corresponded  to  a hypothetical 
case  of  a device  which  detonates  with  its  designed  fission  yield  but 
with  improper  fusion  yield.  The  three  total  yields  considered  were 
2.7  Mt,  3.0  Mt,  and  3.3  Mt,  each  with  a fission  yield  of  1.5  Mt. 
Figure  3. 16  is  the  H+l  hr  normalized  dose  rate  along  the  hotline  for 
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Variations  in  Predicted  Dose  Rate  Due  to  *10  Percent 
Uncertainties  in  Total  Yield 


Uncertainties  in  Total  Yield 
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Figure  3. 16.  Variations  in  Predicted  Dose  Rate  Due  to  *10  Percent 
Uncertainties  in  Fusion  Yield 


the  three  devices,  and  Figure  3. 1?  is  the  total  integrated  hotline 
dose  48  hours  after  burst  for  the  same  devices.  Again,  a maximum 
*10  percent  variation  is  evident. 

Finally,  the  hypothetical  case  of  a device  exploding  with  its 
designed  fusion  yield,  but  not  its  fission  yield,  was  investigated. 

A fusion  yield  of  1.  S Mt  and  fission  yields  of  1. 2 Mt,  1. 5 Mt,  and 
1.8  Mt  were  used  in  this  investigation.  Figure  3. 18  is  a plot  of 
H+l  hr  normalized  dose  rate  for  the  three  devices  along  the  hotline, 
and  Figure  3. 19  is  a plot  of  total  integrated  dose  48  hours  after 
burst  along  the  hotline  for  the  same  devices.  Again,  a maximum 
*10  percent  variation  is  observed. 

These  calculations  have  shown  that  uncertainties  in  the  total 
yield  and  fission  and  fusion  fractions  can  lead  to  uncertainties  of  the 
same  order  of  magnitude  in  the  predicted  dose  rates  and  doses  at 
specific  downwind  ranges. 

3. 4 FISSION  TYPE  SENSITIVITY 

All  the  runs  discussed  in  this  report  were  performed  using 
235 

U high  energy  (U235HE)  as  the  device  fission  type.  The  sensi- 
tivity of  the  results  to  fission  type  are  investigated  in  this  section. 

Tompkins  (Reference  1)  performed  a sensitivity  analysis  of  the 

235 

Particle  Activity  Module  in  which  he  found  that  U fission  spectrum 
239 

(U235FI)  and  Pu  high  energy  (P239HE)  were  the  two  fission  types 
which  bound  the  K-factor  as  a function  of  particle  size  results.  The 
K-factor  relates  exposure  rate  to  debris  deposition  density. 

Two  yields  were  investigated,  30  kt  and  3 Mt.  Figure  3. 20 
is  the  30  kt  dose  rate  along  the  hotline  normalized  to  H+ 1 hr  for  the 
three  fission  types  while  Figure  3.21  is  the  30  kt  total  integrated  dose 
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Variations  in  Predicted  Dose  Rate  Due  to  *10  Percent 
Uncertainties  in  Fission  Yield 
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Figure  3. 20.  Variations  in  Predicted  Dose  Rate  Due  to  Three 
Fission  Types  for  a 30  kt  Device 


Types  for  a 30  kt  Device 


48  hours  after  burst  for  the  same  three  fission  types.  Figures  3. 22 
and  3. 23  are  the  respective  3 Mt  hotline  dose  rates  and  integrated 
doses  for  the  three  fission  types.  It  can  be  seen  that  the  variations 
due  to  fission  type  are  at  most  on  the  order  of  ±10  percent. 

3. 5 SUMMARY 

The  sensitivity  of  fallout  predictions  to  variations  in  total 
yield,  fusion  yield,  fission  yield,  and  fission  type  were  investigated. 
Four  significant  conclusions  were  reached  as  a result  of  this  investi- 
gation. 

First,  variations  in  total  yield,  fission  yield  and  fusion 
yield  are  directly  reflected  in  variations  of  predicted  dose  rate  and 
total  dose  of  the  same  order  of  magnitude  along  the  hotline. 

Second,  changes  in  fission  type  can  affect  dose  rate  and 
dose  predictions  by  no  more  than  ±10  percent. 

Third,  the  maximum  downwind  extent  of  significant  dose 
levels  is  not  a linear  function  of  yield. 

And  last,  to  more  accurately  describe  long-range  fallout, 
many  small  particle  DELFIC  cloud  wafers  should  be  used. 
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Figure  3. 22.  Variations  in  Predicted  Dose  Rate  Due  to  Three 
Fission  Types  for  a 3 Mt  Device 


igure  3. 23,  Variations  in  Predicted  Dose  Rate  Due  to  Three 
Fission  Types  for  a 3 Mt  Device 


Section  4 


SIZE  DISTRIBUTION  SENSITIVITY 


The  subject  of  this  section  is  the  size  distribution  of  the 
entrained  soil  material.  Two  specific  aspects  are  addressed.  The 
first  is  to  show  how  different  analytic  forms  of  the  distribution 
affect  fallout  predictions  for  a surface  burst.  Secondly,  the  height - 
of -burst  dependence  is  discussed. 


SENSITIVITY  TO  ANALYTIC  FORM  OF  THE  SIZE  DISTRI- 
BUTION 


In  the  discussion  to  follow,  the  size  distribution  will  be 
expressed  as  cumulative  mass  fraction  versus  particle  size.  In 
this  form  the  dependent  variable  is  the  fraction  of  the  total  soil  bur- 
den in  the  cloud  with  particles  smaller  than  a specified  diameter. 

The  distribution  refers  to  the  total  amount  of  soil  that  initially  is  in 
the  cloud  and  hence  does  not  include  a time-dependence  due  to  fall- 
out. Section  6 presents  the  results  of  an  investigation  into  the  un- 
certainties associated  with  size-dependence  of  specific  activity 
(defined  as  the  activity  of  a particle  divided  by  its  weight).  Together, 
the  two  distributions  determine  the  relative  sizes  and  positions  and 
shapes  of  iso-dose  contours  on  the  ground. 


Three  size  distributions  were  chosen  for  comparison.  The 
first  two  have  analytic  forms  that  are  power-law,  and  log-normal. 
The  third  is  a hybrid  with  the  low  end  of  the  size  spectrum  having 
log-normal  behavior  while  the  larger  particles  follow  a power-law 


u 


distribution;  the  motivation  for  choosing  these  distributions  is  based 
on  the  following  comments.  Fallout  calculations  using  DELFIC 
and  other  codes  have  traditionally  employed  the  log-normal  distri- 
bution. The  experimental  basis  for  such  a choice  is  found  in  the 
analysis  (Reference  1)  of  cloud  and  fallout  samples  obtained 
at  late  times  (several  hours  after  burst).  Because  of  gravitational 
settling,  the  particles  in  the  samples  tend  to  have  diameters  less 
than  a few  hundred  microns.  The  log-normal  size  distribution  is 


given  by 


where 

CM  = cumulative  mass  fraction 
K = a normalization  constant 

x = log1Q  a 

xo  - lo*10S 

<t  = standard  deviation  of  the  distribution 

a = particle  diameter 

a = median  diameter. 

(2) 

Based  on  an  analysis  of  data  from  shot  SMALL  BOY,  DELFIC 
uses  the  following  values  for  a and  e r unless  the  user  specifies  his 
own  distribution 

a = 0. 407 n 
<r  = logjQ  (4. 0) 

for  this  choice  of  parameters,  50  percent  of  the  mass  is  associated 
with  particles  below  130/*. 
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Complementing  the  studies  of  late -time  cloud  and  fallout 
samples  have  been  investigations  into  the  characteristics  of  soil 
material  ejected  from  the  craters  of  various  tests.  Here  the  particle 
sizes  typically  ranged  from  millimeters  to  meters.  An  analysis^  ^ 
of  the  size  distribution  of  the  crater  ejecta  indicated  a power-law 
behavior,  i.e. , 

CM2(a)  = Ka“P. 

It  was  found  in  reference  3 that  the  various  data  indicated  P to  be 
about  0. 5 for  cohesive  soils  and  rock. 

In  spite  of  the  apparent  disagreement  between  the  log-normal 
and  power -law  distributions,  it  is  important  to  note  that  they  were 
based  on  different  segments  of  the  overall  spectrum  of  particle  sizes 
that  are  initially  in  the  cloud.  In  fact,  a compromise  can  be  reached 
by  creating  a hybrid  distribution,  i.  e. , log-normal  on  the  low  end 
and  power  law  on  the  high  end.  There  are  data  to  support  this 
hypothesis^.  The  analytic  form  of  this  distribution  is 

CMj(a)  = CMj(a)  a < at 

= CMj(at)  + K*  (aP  -aP)  a>at 

where  at  is  diameter  where  the  transition  occurs,  CM^(a)  is  given 
above  and  K'  provides  continuity  at  a^.  Not  explicitly  shown  is  an 
overall  normalization.  For  this  investigation,  the  transition  diam- 
eter was  taken  to  be  upper  limit  on  the  30th  size  class  with  a value 
of  756^i.  Shown  in  Figure  4.1  are  the  three  distributions  as  used 
for  the  300  kt  case.  Note  that  the  log-normal,  the  hybrid,  and  the 
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omparison  of  Size  Distributions 


power-law  cumulative  mass  fractions  attain  50  percent  by  130/i,  460/* 
and  2500/x  respectively. 


As  described  In  previous  sections,  the  base  line  DELFIC 
runs  were  performed  with  a power -law  size  distribution.  The  logic 
of  the  Initial  Conditions  Module  sets  up  the  particle  size  classes  so 
that  each  size  class  contains  an  equal  percentage  of  the  total  mass. 

In  all  the  DELFIC  runs,  40  size  classes  were  chosen,  so  that  the 
mass  fraction  in  each  size  class  was  0. 025.  The  procedure  developed 
to  investigate  the  dependency  of  the  analytic  form  of  the  size  distri- 
bution made  use  of  the  grounded  particles  tape  and  the  logic  of  the 
Initial  Conditions  Module  to  create  new  grounded  particle  tapes  with 
different  size  distribution.  Program  DIAM,  listed  in  Appendix  D, 
was  developed  to  create  a grounded  particles  tape  with  a log-normal 
distribution  from  one  with  a power-law  distribution.  A modification 
of  this  program  was  used  to  create  a tape  with  a hybrid  distribution, 
part  log-normal,  part  power -law. 


The  hybrid  distribution  was  created  with  the  constraint 
that  the  transition  from  log-normal  to  power  law  be  continuous. 
Figure  4. 2 is  a schematic  representation  of  the  process  used.  The 
mass  fractions  generated  by  the  log-normal  distribution  for  the  11 
smallest  size  classes  were  used  in  the  hybrid  distribution  for  all 
yields  along  with  a mass  fraction  of  0.025  for  the  remaining  29 
classes.  These  distributions  were  then  normalized  to  produce  the 
corrected  hybrid  distributions.  Tables  4. 1,  4.2,  and  4.3  are  the 
mass  fractions  in  each  size  class  generated  by  the  power-law,  the 
log-normal,  and  the  hybrid  distributions  for  the  3 Mt,  300  kt,  and 
30  kt  yields,  respectively.  It  can  be  seen  from  these  tables  that 
the  log-normal  distribution  contains  approximately  50  percent  of  the 
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Figure  4.  2.  Schematic  Representation  of  Generation 
of  Hybrid  Distribution 
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Table  4. 1 


3 Mt  Mass  Fractions 


Size  Particle 


lass 

Diameter  (vm) 

Power  Law 

Log-Normal 

Hybrid 

1 

*7*?. 9 1 «* 

.025 

.00014377 

.01552743 

2 

*343. 7**« 

.00014495 

.01552743 

3 

7915. 90' 5 

.00014947 

.01552743 

4 

7499. ,>74  | 

.00021959 

.01552743 

i 

7fl93.9il'-  ) 

.00025251 

.01552743 

* 

**99.  *03'.' 

.00029242 

.01552743 

7 

*31*. 9W 

.00333949 

.0165274  3 

9 

S9<.S.  3 7ji 

.00039519 

.01*52743 

9 

4**6 , ft4-><* 

,0004412* 

.01552743 

10 

S2 45.474? 

.00053997 

.01552743 

11 

4*9*. ( /hS 

.00043397 

.01552743 

12 

4S7| .4*4-> 

,00074445 

.01552743 

13 

4?S*  , 9**<iR 

. 000«5??0 

.01552743 

U 

3952. 32J1 

.001045*9 

.01552743 

IS 

3*S9.SS‘>7 

.00124433 

.01552743 

IS 

3379.051*4 

.00144544 

.01552743 

17 

3107.«073 

.00179119 

.01552743 

1H 

2942.  M2  »*4 

.00214391 

.01552743 

19 

2*01  .ti9-**> 

.00259115 

.0155274  3 

20 

2354  .*  i D 

.00314571 

.01552743 

21 

21  39.4  19* 

.00393**3 

.01659743 

2? 

1925.4*4  < 

.00470197 

,0i55?74J 

23 

172/*. 4 0*2 

.00579145 

.01552743 

2<* 

15  ll . 3«  «.i 

.00717149 

.01692743 

25 

n*i  .22.1* 

.00492974 

.01652743 

2* 

I 1*2. 3 Hi 

.01117957 

.0165274  3 

27 

i ft  24  , * 7** 

.0 1 4079 1 C 

.01552743 

2« 

* 7*».29  1*. 

.017934*1 

.01552743 

29 

743. 1 Ml 

.02272744 

• 0 1 55274  3 

30 

519.304*. 

.02911294 

.01*09190 

31 

50*.7g77 

.03749041 

•0232799b 

32 

405 . 3*42 

.04942*21 

.030077  39 

33 

31S.2W.U 

.06264905 

.03«90**0 

34 

234.4* 1* 

.04076251 

.0501*1  1* 

35 

14*. 4ft  2n 

.10249625 

.043909** 

3* 

l 1 ?. *017 

.12735212 

.07909906 

37 

47.54  U- 

.1475097* 

.09141729 

39 

33. 7 4 j5 

.145425*9 

.09044769 

39 

1 1.2442 

1 

.09145241 

.05*900*6 

40 

1.4075 

.011794** 

.007325*6 
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Table  4.2.  300  kt  Mass  Fractions 


i 


S 


/ 


> 


: 


I 

i 


Size 

Class 

Particle 

Diameter  (jjm)  Power  Law  Log-Normal 

l 

9745. 4i?5  .025  ,00011411 

? 

9?5A.14?0 

•00013129 

i 

8783. 3b64 

.00015139 

<♦ 

81? 1.0830 

.00017499 

5 

7871 ,?947 

.00020275 

b 

7434, OOub 

.0002355? 

7 

7009. ?C0b 

.00027431 

A 

6596.89*7 

•00032035 

9 

6197.0828 

.00037520 

10 

5809.7651 

.00044075 

11 

5434.9416 

.00051934 

12 

5072.61?? 

.00061393 

13 

47??. 7768 

•00072821 

14 

4185.4156 

.00086679 

15 

4060.5885 

.00103556 

16 

3748.2356 

.00124196 

17 

3448.1767 

.0014958? 

1H 

3161.0120 

.00180848 

19 

2886.1414 

.00219664 

20 

2623.7649 

.00268057 

21 

2373.8825 

.00326709 

22 

2136,4943 

.00405154 

23 

1911.600? 

.00502063 

24 

1699.2001 

.00625657 

25 

1499.294? 

.00784263 

26 

1311. 8824 

.00989094 

27 

1136.9648 

•01255319 

28 

974.5413 

.01603518 

29 

824.6118 

.02061645 

30 

687.1765 

.02667523 

31 

5b?.?3S3 

.03471743 

3? 

449.7883 

.0454015? 

33 

349.8363 

.05953363 

34 

2b?.376S 

.07795017 

35 

187.4118 

.10114348 

3b 

124.941? 

.12794149 

37 

74.9647 

.15228020 

38 

37.4824 

.15584745 

39 

12.49*1 

.10303699 

40 

1.5618 

.01430255 

Hybrid 

.01539562 

,01519ms? 

.015395b? 

•015395b? 

.0153956? 

•015395b? 

,015195b? 

,015395b? 

,015195b? 

,015195b? 

,015195b? 

,015395b? 

,015195b? 

,015395b? 

,015195b? 

,015395b? 

,015395b? 

,015195b? 

, 015195b? 

,015395b? 

,015395b? 

.01*3956? 

,015395b? 

,015395b? 

.0151956? 

,015195b? 

,015395b? 

,015195b? 

,015395b? 

.01642727 

• 0? 1 1795b 

.02795939 

,0!bbb??9 

,0480055? 

,0b???bb7 

,0787895b 

.09377794 

.09597974 

,06345?74 

.00880787 


66 


Table  4.3.  30  kt  Mass  Fractions 


Size 

Particle 

Class 

Diameter  (^m) 

Power  Law 

Log-Normal 

Hybrid 

1 

6845.9185 

.025 

.00014153 

.01551792 

2 

8403.0224 

.000162*1 

.01551792 

3 

7972.6675 

.00018679 

.0155179? 

4 

7553. 0o3S 

.00021531 

.01551792 

5 

7J44.7804 

.00024878 

.01551792 

4 

6747.8441 

.00028815 

.01551792 

7 

S.362.2568 

.00033461 

.0155179? 

8 

5988,0064 

.00038959 

.01551792 

0 

5625.0969 

•00045486 

.0155179? 

10 

5273.528* 

.00053258 

.01551792 

11 

4933.3007 

.00062545 

.0155179? 

12 

4604.4140 

.00073680 

.0155179? 

13 

4286.868? 

.00087080 

.01551792 

1 4 

3980  • 64.14 

.00103266 

.01551792 

IS 

3685.79-44 

.00122694 

.01551792 

14 

3402.2744 

.00146795 

.P1S5179? 

17 

3130.094? 

.00176023 

.01551792 

IS 

?«69.?5  11 

.00211927 

.01551792 

19 

2619.7526 

.00256236 

.0155179? 

20 

2381 .59  14 

.00311186 

.0155179? 

21 

2154.7750 

.00379676 

.01551792 

22 

1939.2975 

. 00*65494 

.0155179? 

23 

1735.1609 

.00573609 

.0155179? 

24 

1542. 36o3 

.00710581 

.0155179? 

25 

1360.9105 

.00885110 

.01551792 

26 

1 1 90.7947 

.01108785 

.01551792 

27 

1032.0238 

.01397104 

.01551792 

28 

884 ,5919 

.01770808 

.01551792 

29 

748,5008 

.02257624 

•C 155 1 79? 

30 

623.7507 

.02894357 

.01796575 

31 

510.3414 

.03729033 

.0231*672 

32 

408 .2712 

.04822043 

.029931 35 

33 

317 .5456 

.06243442 

.03875409 

34 

238.1893 

.08058091 

.05001791 

35 

170.1 138 

.10279304 

.06380535 

36 

1 13.4092 

.12741502 

.0790886? 

37 

68 . 0465 

.14785776 

.09177777 

38 

34.0228 

.14633127 

.09081025 

39 

1 1.3409 

' 

.09222029 

. 05724266 

40 

1.4176 

.01195390 

.00741998 

67 


mass  in  the  five  smallest  size  classes,  and  that  the  normalization 
of  the  hybrid  distribution  results  in  there  being  approximately  60 
percent  of  the  mass  in  each  class  with  respect  to  the  "parent" 
distributions.  These  changes  in  the  relative  mass  distributions 
directly  affect  the  dose  rates  predicted  on  the  ground. 

Figure  4. 3 presents  the  3 Mt,  H+l  hr  normalized  dose  rate 
for  the  three  distributions  along  a line  345°  from  the  horizontal  (the 
chosen  "hotline")  as  a function  of  range.  The  interpretation  of  these 
curves  is  quite  simple  in  terms  of  aerodynamics.  The  smaller  par- 
ticles will  be  raised  higher  by  the  buoyant  cloud,  will  fall  slower 
through  a viscous  atmosphere,  and  consequently,  will  impact  the 
ground  later  and  farther  along  the  downwind  axis  than  the  larger 
particles.  Impact  time  and  distance  along  the  downwind  axis  are 
directly  related  to  particle  diameter. 

The  higher  dose  rates  predicted  in  the  vicinity  of  ground 
zero  by  the  power -law  distribution  are  caused  because  of  the  greater 
percentage  of  mass  allocated  to  the  larger  particles.  That  this  dose 
rate  is  directly  proportional  to  this  mass  fraction  can  be  seen  from 
a comparison  of  the  fractions  in  the  large  particle  sizes  for  the 
power-law  and  hybrid  distributions.  The  power-law  fractions  are 
almost  one  and  a half  times  the  hybrid  fractions,  and  this  is  the 
ratio  of  the  two  dose  rates  for  approximately  the  first  20  km  down- 
wind. In  the  same  way,  it  can  be  readily  observed  that  the  par- 
ticles impacting  from  approximately  30  km  downwind  and  further 
are  in  the  11  smallest  size  classes  since  the  hybrid  curve  values 
are  two-thirds  of  the  log-normal  values  in  this  region. 


These  curves  show  that  variations  in  the  size  (and  therefore  ' 

mass)  distributions  can  cause  variations  in  the  predicted  dose  on  the  j 

order  of  a factor  of  20  near  ground  zero,  and  a factor  of  5 at  larger  i 

ranges.  Data  are  only  shown  along  the  hotline  to  140  km  because 
past  this  range  the  probability  of  single  wafer  impact  is  greater,  and 
die  effects  will  then  be  proportional  to  mass  fractions  in  individual 
size  classes  and  not  to  overall  size  distributions. 

Figure  4.4  presents  the  300  kt,  H+l  hr  normalized  dose 
rate  for  the  three  distributions  along  a line  345°  from  the  horizontal 
as  a function  of  range.  The  same  conclusions  that  were  drawn  from 
the  3 Mt  curves  can  also  be  drawn  here,  the  only  difference  being 
that  the  larger  particles  do  not  travel  as  far  as  in  the  3 Mt  case. 

This  was  the  only  yield  for  which  the  horizontal  wafer  subdivision 
option  was  used,  which  explains  the  relative  roughness  of  this  curve. 

However,  the  factor  of  20  difference  near  ground  zero  and  the  factor 
of  5 difference  further  downwind  are  still  clearly  defined  for  the 
power-law  and  log-normal  distributions. 

Figure  4.  5 presents  the  30  kt,  H+l  hr  normalized  dose 
rate  for  the  three  distributions  along  a line  350°  from  the  horizontal 
as  a function  of  range.  Again,  the  same  effects  that  were  seen  for 
the  larger  yields  are  also  shown  here. 

This  investigation  of  the  fallout  sensitivity  to  the  analytic 
form  of  the  size  distribution  of  dust  has  shown  that  the  overall 
variation  of  the  predicted  dose  rate  can  be  as  large  as  a factor  of  20 
near  ground  zero  and  as  large  as  a factor  of  5 further  downwind  for 
two  reasonable  forms  of  the  size  distribution.  This  occurs  because 
the  predicted  dose  is  proportional  to  the  relative  amounts  of  mass  in 
each  size  class,  and  the  particle  impact  positions  are  dependent  on 
particle  diameter,  not  total  amount  of  mass  at  that  diameter. 
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4. 2 HEIGHT -OF-BURST  EFFECTS  ON  FALLOUT 

This  section  deals  with  the  height -of -burst  (HOB)  dependence 
on  that  part  of  the  fallout  hazard  that  is  attributed  to  fission  debris. 
The  variation  with  HOB  of  the  neutron -activated  soil  is  discussed  in 
Section  7. 

There  are  presently  two  algorithms  for  predicting  the 
change  in  the  fallout  pattern  for  a low -air  burst  versus  a surface 
burst.  Both  determine  a correction  factor  versus  HOB  that  is  to 
be  applied  uniformly  over  the  entire  fallout  pattern  of  a surface 
burst.  The  first  correction  factor  is  a fit  to  empirical  data.  For 
a number  of  U.  S.  tests  at  varying  heights  of  burst,  an  integration 
was  performed  over  the  fallout  pattern  to  determine  the  total  activity 
brought  down  in  the  so-called  "close-in  fallout"  (which  usually 
Implies  times  of  arrival  up  to  H+24  hrs)  . The  form  of  this  correc- 
tion factor^  is 


Ft  = exp  (- .01238  X) 

where 

X , is  the  scaled  HOB,  given  by 


X = HOB  -W 

with  the  HOB  in  feet  and  W,  the  yeild,  in  kt.  There  is  a large  scatter 
in  the  data  on  which  the  fit  is  based  so  that  the  uncertainty  varies 
from  *50  percent  at  A = 50  feet  to* an  order  of  magnitude  atx=  300  feet. 

The  second  correction  factor^  is  based  on  the  volume  of 
intersection  between  a sphere  and  the  ground.  The  radius  of  the 
sphere  is  given  by 
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R = 180WU>* 

where  R is  in  feet.  The  form  of  the  correction  factor  is 

F2  = 0.  5 (1  - HOB/R)2  (2  + HOB/R). 

Models  such  as  these  omit  potentially  significant  effects  that  change 
die  size  distribution  of  the  radioactive  soil  and  thus  result  in  a fall- 
out pattern  that  cannot  be  related  to  that  of  a surface  burst  by  a sim- 
ple correction  factor. 

In  the  case  of  a surface  burst,  most  of  the  dust  is  crater 
material  ejected  into  the  fireball  before  rise  begins  and  there  is 
probably  homogeneous  mixing  between  debris  vapor  and  dust,  part 
of  which  may  itself  be  melted  or  vaporized.  (Further  research 
would  be  required  to  confirm  this  observation,  however. ) As  the 
HOB  increases,  however,  two  changes  occur  which  make  the  estima- 
tion of  fallout  dose  levels  more  difficult.  First,  at  a low  HOB 
(approximately  15  feet  for  a 1 kt  yield),  the  ejection  crater  disappears 
although  compaction  might  still  occur.  At  HOBs  above  the  crater - 
formation  cutoff,  the  one  remaining  source  of  dust  is  the  thermally 
generated  surface  layer.  Entrainment  of  this  material  is  not  instan- 
taneous. Depending  on  the  HOB,  it  may  occur  after  the  debris  has 
cooled  and  condensed  onto  whatever  nucleation  sites  are  normally 
available  in  the  atmosphere.  Furthermore,  the  spatial  distribution 
of  the  debris  is  changing  from  symmetrical  to  toroidal. 

The  approach  taken  to  examine  HOB  effects  employed  the 
VORDUM  code  described  elsewhere  in  this  report  to  determine  the 
fraction  of  the  debris  mass  that  was  co-located  with  dust  of  various 
sizes  at  several  times.  The  source  of  the  dust  is  the  surface  layer. 
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The  bomb  debris  was  modeled  as  flow  field  trace  particles  distributed 
in  a uniform  density  cylinder  with  radius  equal  to  0. 5 the  fireball 
radius.  The  vertical  cylinder  was  fixed  so  that  its  volume  equaled 
that  of  a sphere  with  a radius  half  the  fireball  radius.  These  assump- 
tions are  based  on  the  best  current  estimation  for  describing  the 
bomb  debris.  There  is,  however,  continuing  controversy  about  the 
radial  extent  of  the  debris  at  the  time  of  pressure  equilibrium.  One- 
dimensional radiation  hydrocodes  predict  the  debris  to  be  within  0. 1 
fireball  radii  of  the  burst  point  while  photographic  evidence  indicates 
0. 5 radii.  Currently  it  is  the  feeling  of  some  researchers  in  this 
area  that  the  hydrocode  calculations  are  incomplete  in  that  turbulence 
is  not  included.  Figure  4. 6 is  a schematic  of  the  two  sources  for  a 
typical  configuration. 

Because  the  time  at  which  fission  debris  condenses  varies 
depending  on  the  volatility  of  each  mass  chain,  a range  of  times 
was  selected  for  examination  — 5 to  30  seconds.  No  attempt  is 
made  in  this  study  to  model  the  temperature  history  of  the  debris. 

In  the  VORDUM  calculation,  the  debris  and  the  dust  are 
represented  as  sets  of  trace  particles.  Associated  with  each  trace 
particle  is  a fraction  of  the  total  debris  or  dust  mass.  The  dust 
trace  particles  are  tracked  through  the  flow  field  and  at  several 
times  a snapshot  of  their  positions  is  obtained  which  defines  a two- 
dimensional  envelope  of  the  dust  region.  The  debris  trace  particles 
are  tracked  through  the  same  flow  field.  At  each  snapshot  time. 


their  individual  locations  are  examined  to  determine  which  (if  any) 
are  within  the  dust  region.  If  a particular  debris  trace  particle  is 
within  the  dust,  then  the  fraction  of  total  debris  which  it  represents 
is  defined  to  be  mixed  with  the  dust.  Shown  in  Figure  4. 7 is  a 
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Figure  4. 6.  Schematic  of  Dust  and  Debris  Sources 
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Figure  4.  7.  Comparison  of  Dust  and  Debris  Regions 
for  a Low  Air  Burst 
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typical  output  from  VORDUM  using  the  two  sources  described  above. 
The  parameters  were  a 300  kt  burst  at  1050  feet  and  the  time  at  15 
seconds.  The  solid  contour  outlines  the  debris  region  while  the  two 
dashed  contours  correspond  to  dust  particles  with  diameter  of  0.025 
cm  and  0. 1 cm  respectively.  The  significance  of  this  figure  is  that 
if  15  seconds  is  the  time  at  which  the  debris  condenses,  then  a pre- 
diction can  be  made  about  the  fraction  of  the  total  debris  that  is  not 
co-located  with  dust  particles  of  a particular  size.  In  this  way,  an 
upper  limit  can  be  obtained  for  the  fraction  of  the  debris  that  con- 
denses onto  dust  particles  as  a function  of  size  and  HOB.  Determi- 
nation of  such  an  upper  limit  avoids  difficulties  related  to  the  conden- 
sation process  itself  as  well  as  fractionation  over  the  size  spec- 
trum of  the  dust.  Hence,  this  analysis  provides  only  a partial 
answer  to  the  determination  of  HOB  effects.  As  will  be  seen  below, 
significant  differences  exist  between  these  upper  limits  and  the 
algorithms  presently  used  in  fallout  codes  and  discussed  above. 

Three  yields  were  investigated,  50  kt,  200  kt  and  1 Mt.  The 
VORDUM  code  was  run  for  each  at  several  heights  of  burst  ranging 
from  96  feet  (scaled  to  1 kt  by  the  cube  root  of  the  yield)  up  to  275 
feet.  No  runs  were  made  for  cratering  bursts.  Four  particle  sizes 
were  used  — 0.025  cm,  0.050  cm,  0. 1 cm,  and  0.  5 cm.  The  small- 
est of  these  is  run  in  VORDUM  with  the  assumptions  of  infinite  drag 
coefficient  and  zero  gravity  so  that  the  flow  field  is  tracked  perfectly. 
Consequently  the  0.025  cm  trace  particles  represent  particles 
smaller  than  this  size.  Figures  4.8  - 4. 10  show  the  HOB  (scaled 
to  1 kt)  dependence,  of  the  fraction  of  debris  mixed  with  dust 
(up  to  0. 025  cm)  for  the  three  yields  respectively.  Also,  indicated 
on  each  figure  are  the  two  algorithms  — F^  is  a fit  to  test  data  and 
F2  is  based  on  sphere-ground  intersection.  In  the  50  kt  case,  it  is 


Fraction  Relative  to  Surface  Burst 


1.0 


Figure  4.8.  Fraction  of  Debris  Mixed  with  Dust  (Less  Than  .025  cm) 
Versus  HOB  for  50  kt 
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Figure  4. 9.  Fraction  of  Debris  Mixed  with  Dust  (Less  Than 
.025  cm)  Versus  HOB  for  200  kt 
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Fraction  Relative  to  Surface  Burst 
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W = 1 MT 

Dust  Size  < .025  cm 
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Figure  4. 10.  Fraction  of  Debris  Mixed  with  Dust  (Less  loan  .025  cm) 
Versus  HOB  for  1 Mt 
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apparent  that  is  overestimating  the  fallout  activity  for  scaled 
HOBs  above  200  feet.  The  present  calculations  tend  to  support  the 
Fg  curve  for  times  greater  than  15  seconds.  However,  if  debris 
condensation  and  attachment  to  the  dust  is  occurring  much  earlier 
(e.  g. , 5 sec),  then  Fg  is  also  an  overestimate.  Comparison  at 
200  kt  again  indicates  F1  to  be  higher  than  either  F2  or  the 
calculations  above  a scaled  HOB  of  275  feet.  At  1 Mt  a similar  com- 
parison can  be  made. 

Figures  4. 11  through  4. 18  depict  Fj,  Fg  and  the  calculations 
of  dust-debris  mixing  for  the  three  yields,  but  for  dust  particle  diam- 
eters from  0. 05  cm  to  0. 5 cm.  The  dust  trace  particles  used  for 
these  sizes  were  subject  to  finite  drag  and  gravity.  The  fraction  of 
debris  mixed  with  dust  decreases  with  increasing  size.  For  example 
Figures  4. 11  and  4. 14  show,  in  contrast  with  the  small  dust  sizes 
(Figure  4.  8)  that  for  50  kt,  both  the  Fj  and  F2  algorithms  are  too 
high  for  scaled  HOBs  above  100  feet.  However,  this  size  effect 
has  a yield  dependence  as  can  be  seen  from  Figures  4. 13,  4. 16  and 
4.18. 

An  important  observation  on  the  calculations  is  the  HOB  above 
which  there  is  no  mixing  by  30  seconds.  Although  not  shown  on  the 
figures,  all  yields  were  examined  for  the  HOB  such  that  the  fireball 
just  touches  the  ground.  In  no  case  was  there  any  mixing  of  dust 
( of  any  size)  and  debris. 

The  above  analysis  is  incomplete  in  the  sense  that  the 
condensation  of  the  fission  debris  is  not  modeled.  However,  the 
results  do  indicate  that  for  yields  above  200  kt,  the  Fg  algorithm 
is  more  realistic  than  F^  for  scaled  heights  of  burst  above  200  feet. 
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Figure  4. 12.  Fraction  of  Debris  Mixed  with  Dust  (0. 5 cm) 
Versus  HOB  for  1 Mt 
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Figure  4. 13.  Fraction  of  Debris  Mixed  with  Dust  (0.05  cm) 
Versus  HOB  for  50  kt 
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Figure  4.14.  Fraction  of  Debris  Mixed  with  Dust  (0.05  cm) 
Versus  HOB  for  200  kt 


Figure  4. 15.  Fraction  of  Debris  Mixed  with  Dust  (.05  cm) 
Versus  HOB  for  1 Mt 
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Figure  4.  la  Fraction  of  Debris  Mixed  with  Dust  (0. 1 cm) 
Versus  HOB  for  50  kt 
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Both  of  these  prescriptions  overestimate  by  approximately  an  order 
of  magnitude.  The  very  close-in  fallout  (due  to  dust  larger  than 

0.05  cm  which  will  have  times  of  arrival  up  to  an  hour)  from  low 
yield  (e.g. , 50  kt)  weapons  burst  above  150  scaled  feet. 


4.3  REFERENCES 

1.  R.  Heft,  "The  Characterization  of  Radioactive  Particles 
from  Nuclear  Weapons  Tests,"  UCRL-70790,  Lawrence 
Livermore  Laboratory,  University  of  California,  Liver- 
more, California,  1968. 

2.  J.  Maloney,  Ballistic  Research  Laboratories,  U.  S.  Army, 
Aberdeen  Research  and  Development  Center,  Aberdeen, 
Maryland,  private  communication,  1974. 

3.  W.  M.  Layson,  Science  Applications,  Inc. , McLean,  Vir- 
ginia, Private  communications,  1974. 

4.  M.  Tobriner,  et  al. , "Technology  Support  for  Environmental 
Defense  — Second  Bimonthly  Progress  Report, " SAI-74- 
561-WA,  Science  Applications,  Inc.,  McLean,  Virginia,  1974. 

5.  R.  L.  Showers,  "Improvements  to  the  PROFET  Fallout 
Prediction  Program, " BRL  Memorandum  Report  2095, 
(AD-883-280),  Ballistic  Research  Laboratories,  U.  S. 

Army  Aberdeen  Research  and  Development  Center, 

Aberdeen,  Maryland,  February,  1971. 

6.  Draft,  Effects  of  Nuclear  Weapons,  October,  1973. 


91 


Section  5 


SENSITIVITY  TO  CLOUD  MODEL  ASSUMPTIONS  DURING  RISE 

5. 1 NONUNIFORM  MASS  DISTRIBUTION 

The  sensitivity  of  the  fallout  predictions  to  the  DELFIC 
assumption  of  uniform  density  mass  wafers  was  investigated  using 
the  output  from  a 3 Mt  VORDUM  run.  VORDUM  is  an  analytic 
model  which  has  been  developed  to  describe  the  air  and  dust/debris 
motion  inside  and  around  a nuclear  cloud.  Inside  the  cloud  the  air 
motion  takes  the  form  of  a spheroidal  vortex  which  persists  as  the 
cloud  rises,  while  the  motion  outside  the  cloud  is  that  induced  by  a 
spheroid  moving  through  an  inviscid  fluid.  The  model  uses  empiri- 
cal equations  for  rise  velocity  and  radial  expansion. 

The  procedure  followed  in  making  comparisons  was  to  first 
insert  the  DELFIC  rise  and  expansion  equations  into  VORDUM.  A 
total  of  1656  trace  particles,  representing  the  same  particle  size 
class  range  that  was  used  for  DELFIC,  was  run  in  VORDUM  until 
stabilization  time.  The  cloud  defined  by  these  particles  was  then 
radially  expanded  according  to  the  DELFIC  equations  until  848  seconds 
after  burst,  the  time  at  which  the  Cloud  Rise  Module  ends.  Figure  5. 1 
is  a comparison  of  the  two  clouds  at  this  time.  Since  both  VORDUM 
and  DELFIC  produce  cylindrically  symmetric  clouds,  this  figure 
consists  of  cloud  contours  in  the  R and  Z coordinate  system.  Shown 
are  the  maximum  cloud  extents  and  the  average  densities  within  the 
clouds.  The  DELFIC  cloud  is  larger  radially  and  contains  less  mass 
than  the  VORDUM  cloud,  and  is  consequently  less  dense.  The 
VORDUM  cloud  also  contains  larger  particles. 
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A comment  on  the  generation  of  the  above  cloud  r ..  tours  is  in 
order.  VORDUM  has  associated  with  it  a density  program  which 
outputs  mass  and  number  densities,  as  well  as  radar  reflectivity  and 
attenuation  information.  No  such  program  exists  with  DELFIC. 
Subroutine  RSXP  has  the  option  of  printing  the  pertinent  information 
for  each  wafer,  and  this  information  was  used  in  producing  the 
desired  contours.  However,  an  examination  of  the  output  showed 
certain  discrepancies,  and  modifications  had  to  be  made  to  produce 
meaningful  results.  Appendix  A contains  an  explanation  of  these 
discrepancies  and  the  methods  used  to  correct  them. 

Comparisons  between  the  particle  positions  predicted  by  the 
two  codes  were  made  using  the  procedures  discussed  above.  Figures 
5. 2 through  5.  5 show  the  volumes  of  space  occupied  by  the  specified 
particles  transported  through  the  two  codes.  It  can  be  seen  that 
overlapping  volumes  exist  only  for  particles  smaller  than  0. 1441  cm. 

Discrepancies  in  the  grounded  particles  tape,  discussed  in 
Appendix  B,  as  well  as  the  large  non -overlapping  regions  of  the  two 
clouds,  make  a meaningful  comparison  of  the  fallout  from  the  two 
models  extremely  difficult.  The  following  simplifying  assumptions 
were  made  to  reduce  the  degree  of  difficulty: 

a.  Only  wafers  aloft  at  848  seconds  were  considered. 

b.  Only  wafers  co-located  with  similar  VORDUM 
particles  were  considered. 

c.  The  mass  in  each  wafer  was  constant,  but  the  den- 
sity was  changed  to  reflect  tne  generally  smaller 
VORDUM  particle  occupied  regions. 
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Figure  5.2.  Comparison  of  Particle  Positions 
Produced  by  VORDUM  and  DELFIC 
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Figure  5.3.  Comparison  of  Particle  Positions 
Produced  by  VORDUM  and  DELFIC 
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d.  Only  VORDUM  particle  positions  were  important, 
not  relative  density. 

e.  The  modified  DELFIC  wafer  mass  was  still  uni- 
formly distributed  over  the  annuli. 

Figure  5. 6 is  a simplistic  example  of  the  procedure  used  in 
generating  a new  grounded  particle  tape.  The  3 Mt  DELFIC  tape, 
which  contained  1987  particles,  was  modified  to  one  of  323  particles. 
These  particles  represent  most  of  the  long  range  fallout  particles. 

A modified  version  of  LINTAP  (discussed  in  Section  2)  was  produced 
to  predict  fallout  from  the  annuli  created  from  the  wafers  so  that 
one-to-one  comparisons  could  be  made. 

Figures  5.  7 through  5. 9 are  dose  rate  normalized  to  H+l 
hour,  integrated  dose,  and  time  of  arrival,  respectively,  of  the 
fallout  predicted  from  both  clouds  along  a line  345°  to  the  horizontal, 
the  chosen  "hotline. " These  curves  can  be  interpreted  quite  easily 
by  reference  to  Figure  5. 1.  The  radius  of  the  VORDUM  cloud  is 
approximately  two-thirds  that  of  the  DELFIC  cloud,  so  if  the  den- 
sities in  each  wafer  were  changed  to  reflect  just  this  fact,  the  dose 
rate,  which  is  proportional  to  the  density,  should  be  higher  by  a 
factor  of  2 for  the  VORDUM  density  cloud.  Also,  since  the  VORDUM 
modified  wafers  are  generally  smaller  than  the  standard  wafers, 
they  do  not  extend  as  far  along  the  hotline,  and  consequently,  the 
time  of  arrival  more  closely  reflects  the  wafers  which  fell  earlier 
uprange  and  extend  over  a wide  area. 

In  summary,  it  appears  that  the  assumptions  of  uniform 
density  mass  wafers  can  vary  the  predicted  fatlout  by  a factor  of  2, 
and  the  time  of  arrival  calculations  can  be  1000  seconds  early. 
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1)  Volumes  Occupied  by  VORDUM  and  DELFIC  Particles  of 
Diameter  d,  at  848  seconds  after  burst. 

3)  Only  Wafers  1(  2,  3,  and  4 retained. 

3)  Total  Mass  Aloft  = Mass#1  + Mass^  + Mass#j  + Mass^ 


4) 


5) 


New  Wafer  Radii 
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Figure  5. 6.  Pictorial  Representation  of  Simplifying  Assumptions 
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Figure  5. 8.  Comparison  Between  Uniform  Density  Wafers 
and  VORDUM  Mass  Distribution  Annuli  Along 
the  Cloud  Hotline 
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Figure  5. 9.  Comparison  Between  Uniform  Density  Wafers 
and  VORDUM  Mass  Distribution  Annuli  Along 
the  Cloud  Hotline 
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SENSITIVITY  TO  CLOUD  RISE  DYNAMICS 


DELFIC  calculates  the  rise  history  of  the  cloud  thermo- 
dynanxically  including  the  effects,  for  example,  of  latent  heat  release 
as  the  entrained  atmospheric  water  vapor  in  the  cloud  condenses. 

The  principal  drawback  of  this  approach  is  that  some  parameters 
must  be  predetermined  since  there  is  at  least  one  more  unknown  than 
constraining  equations.  In  particular,  the  cloud's  horizontal  radius 
is  obtained  from  the  volume  and  vertical  thickness  assuming  an 
oblate  spheroidal  shape.  VORDUM,  on  the  other  hand,  uses  empiri- 
cally determined  equations  for  rise  velocity,  radius,  and  cloud  top. 

At  high  yields,  the  standard  deviation  on  fits  to  these  parameters  is 
typically  on  the  order  of  10  percent.  In  the  low-kiloton  range,  how- 
ever, there  is  significantly  higher  uncertainty  since  atmospheric 
conditions  especially  the  lapse  rate  play  a more  significant  role. 

In  the  course  of  this  study,  it  was  found  that  the  DELFIC 
and  VORDUM  rise  parameters  differed  significantly.  To  determine 
the  effect  this  has  on  the  fallout  calculations,  VORDUM  was  exercised 
for  three  yields  - 30  kt,  300  kt,  and  3 Mt,  using  both  the  standard 
VORDUM  equations  for  rise  and  expansion  and  the  DELFIC  predictions 
for  these  parameters.  Figure  5. 10  is  a comparison  of  cloud  shapes 
and  densities  for  the  standard  VORDUM  and  the  special  case  with 
DELFIC  rise  parameters  at  4.75  minutes  for  the  300  kt  case.  It 
can  be  seen  that  the  standard  VORDUM  lofts  the  dust  higher  but  the 
two  clouds  are  similar  in  shape  and  structure. 

To  make  a comparison  on  the  basis  of  fallout  dose  levels, 
the  VORDUM  program  was  modified  to  include  the  same  wind  field 
used  in  DELFIC  for  this  study.  In  addition,  VORDUM  was  extended 
in  time  to  calculate  the  time  and  position  of  grounding  for  the  trace 
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Figure  5. 10.  Comparison  of  Modified  VORDUM  Cloud  with  Standa: 
Cloud  for  W = 300  kt,  at  T = 4. 75  Minutes 


particles.  A modified  version  of  LINTAP  was  created  to  calculate 
the  predicted  dose  rates  for  the  two  sets  of  VORDUM  trace  particles 
at  each  of  the  three  yields.  The  major  modification  was  the  assign- 
ment of  mass  to  annuli  rather  than  to  disks.  A grounded  particles 
tape  was  created  for  input  to  this  modified  LINTAP. 

Figure  5. 11  compares  the  standard  VORDUM  rise  with  that 
predicted  by  the  Cloud  Rise  Module  for  30  kt.  The  cloud  bottom  pre- 
dicted by  VORDUM  stabilizes  later  and  much  higher  than  DELFIC 
would  predict  (8.  7 km  versus  5 km).  Figure  5. 12  shows  the  dose 
rate  (normalized  to  H+l  hr)  that  results  using  the  two  rise  formalisms 
in  VORDUM.  The  irregularity  in  the  curves  is  due  to  the  limited 
number  c:  trace  particles  used  in  the  VORDUM  runs  but  it  would 
appear  that  the  average  deviation  for  the  range  0 to  8 km  is  on  the 
order  of  a factor  of  5 to  10.  However,  local  variations  can  be  several 
orders  of  magnitude. 

At  higher  yields,  the  discrepancy  in  rise  is  not  quite  so 
severe.  Figure  5. 13  depicts  the  cloud  bottom  and  radius  for  the 
300  kt  case  and  Figure  5. 14  compares  the  dose  rate  at  H+l  hr  for 
the  standard  and  the  DELFIC -modified  cases.  As  can  be  seen,  there 
is  significantly  better  agreement  here  than  at  lower  yields.  Figure 
5. 15  shows  the  comparison  for  3 Mt  and  confirms  this  observation. 

5.3  REFERENCES 

1.  J.  Shannon,  Science  Applications,  Inc. , McLean, 
Virginia,  Private  communication,  1974. 
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30  kt 


Figure  5.11.  Comparison  of  Rise  Parameters  for  30  kt 


Dose  Rate  (R/Hr) 
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Figure  5. 12,  Comparisons  of  Predicted  Dose  Rate  on  Cloud 
"Hotline”  Normalized  to  H+l  Hour  for  a 30  kt 
Burst 
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Figure  5. 13.  Comparison  of  Rise  Parameters  for  300  kt 


Figure  5. 14.  Comparisons  of  Predicted  Dose  Rate  on  Cloud 
"Hotline"  Normalized  to  H+l  Hour  for  a 300  kt 
Burst 
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Figure  5. 15.  Comparisons  of  Predicted  Dose  Rate  on  Cloud 
"Hotline"  Normalized  to  H+l  Hour  for  a 
3 Mt  Burst 


\ 

C' 

v 


111 


Section  6 


SENSITIVITY  TO  ASSUMPTIONS 
IN  THE  PARTICLE  ACTIVITY  MODULE 


The  subject  of  this  section  is  the  assignment  of  fission  activ- 
ity to  the  dust  particles.  The  Particle  Activity  Module*  uses  the 
radial-distribution  model  of  Freiling  modified  to  account  for  an 
observed  trend  to  constant  specific  activity  for  particles  above  a 
few  hundred  microns. 

For  a given  mass  chain,  the  Freiling  ratio,  FR,  is  defined 
as  the  fraction  of  its  membership  that  has  a higher  melting  tempera- 
ture than  the  entrained  soil  and  hence  is  considered  refractory.  The 
specific  activity  (activity  per  unit  mass)  is  then  assumed  to  be  pro- 
portional to  the  (b.  - 1)  power  of  the  particle  size  where  b.  is  the 

* th  * 

square  root  of  FR  for  the  i mass  chain.  The  reader  is  referred  to 

Reference  1 for  a comprehensive  discussion  of  the  Freiling  model. 


6. 1 ALTERNATIVE  FORM  OF  THE  FREILING  MODEL 

A more  exact  formalism  of  the  Freiling  model  was  examined 
for  its  effect.  In  this  alternative  approach,  the  refractory  and  volatile 
elements  of  each  mass  chain  are  separate.  The  number  of  equivalent 
fissions  of  mass  chain  i in  size  class  k is  given  by 


F.  (dk)  = FT  • Y.  • f (dk) 
i l m 


FR  + 


(1  - FR)  dk 


-1 


. (dk)  dk 
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where  FT  is  the  total  number  of  fissions  in  all  size  classes,  Y.  is 

th  * 

the  fission  yield  of  the  i mass  chain,  fm(dk)  the  fraction  of 

fissions  in  the  k**1  size  class  for  a perfectly  refractory  mass  chain 

and  dk  is  the  mean  diameter  of  the  k**1  size  class.  This  equation. 

however,  has  to  be  modified  to  account  for  the  experimental  evidence 

that  Fj(dk)  becomes  constant  for  large  particles.  Following  the  same 

procedure  described  in  Reference  1,  then. 


Fj  (dk)  = FT  Yj  f m (dk) 


(dk)  dk 


[l.O-Rjd  - 

HI 


where  is  the  fraction  of  fissions  in  the  i mass  chain  that  follows 
the  radial  distribution;  it  is  determined  as  in  the  present  version  by 
that  particle  diameter,  D,  where  the  radial  and  the  constant  specif- 
ic activity  components  cross  (see,  for  example,  Figure  2. 1 of 
Reference  1).  The  fraction,  R. , is  thus  determined  to  be 


Ri  ■ (*• 


0 + FR  + (1  - FR) 


-)"• 


Currently  D is  100m  in  the  PAM.  Shown  in  Figure  6.1  is  a compari- 
son of  the  activity  distribution  with  size  for  the  current  formalism 
(labeled  PAM)  and  that  described  above  for  a 3 Mt  surface  burst. 

The  significance  of  the  comparison  is  that  the  exact  form  predicts 
the  activity  associated  with  particles  greater  than  a few  hundred  mi- 
crons to  be  about  50  percent  lower  than  the  current  formalism  while 
the  activity  on  the  very  small  particles  (less  than  10m)  is  proportionally 
higher.  Figures  6. 2 and  6. 3 show  the  integrated  dose  along  the  hotline  for 
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Comparison  of  Activity  Versus  Size  for  Alternative 
Forms  of  Freiling  Model 
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Figure  6. 2.  Comparison  of  Alternative  Freiling  Model  Forms 
for  3 Mt  Hotline  Dose  - Short  Range 
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Figure  8. 3.  Comparison  of  Alternative  Freiling  Model  Forms 
for  3 Mt  Hotline  Dose  - Intermediate  Range 
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this  yield  in  two  range  regimes,  the  first  extending  to  28  km  and 
the  second  to  280  km.  The  dose  predicted  by  the  standard  formalism 
is  always  larger  typically  by  30  to  50  percent,  but  the  two  become 
comparable  past  160  km.  Not  shown,  are  the  long  range  fallout, 
where,  as  indicated  by  the  comparison  shown  on  Figure  6. 1,  the 
alternative  form  of  the  Freiling  model  would  predict  higher  dose 
levels.  However,  at  these  distances,  the  dose  is  small  enough  that 
50  percent  effects  are  unimportant. 

A comparable  situation  is  manifest  at  other  yields.  As 
shown  in  Figures  6. 4 and  6. 5 for  30  kt,  the  alternative  form  of  the 
Freiling  model  results  in  a lower  dose  out  to  70  km,  but  only  between  5 
and  7 km  is  the  difference  more  than  a factor  of  50  percent.  The 
conclusion  to  be  reached  is  that  while  the  present  coding  is  not  the 
most  exact,  it  overestimates  the  integrated  dose  by  factors  less 
than  50  percent  for  particles  above  20-30  microns  (i.e.,  for  close- 
in  fallout). 

In  addition  to  establishing  the  sensitivity  to  alternative  forms 
of  the  Freiling  model,  the  effect  of  variations  in  the  crossover 
diameter,  D,  was  also  investigated.  This  is  the  diameter  at  which 
the  radially  dependent  and  constant  components  of  the  specific 
activity  cross.  Presently,  it  is  set  at  lOOn.  Shown  in  Figures  6.6 
and  6. 7 are  the  standard  curve  for  integrated  dose  along  the  hot- 
line for  3 Mt  as  well  as  the  curves  that  result  when  D is  set  to 
20u  and  500^ . As  is  indicated  on  these  figures,  there  is  very  little 
sensitivity  to  variations  in  this  parameter  with  the  two  extreme 
values  of  D resulting  in  an  average  variation  on  the  dose  of  about 
plus  or  minus  10  percent.  Up  to  a range  of  200  bn,  the  20 u curve 
is  higher  than  the  standard  and  the  500*/  curve,  bat  past  260  km,  it 
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Figure  6. 4.  Comparison  of  Alternative  Freiling  Model  Forms 
for  30  kt  Hotline  Dose  - Short  Range 
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{ Figure  6. 6.  Sensitivity  of  Hotline  Integrated  Dose  to  Crossover 

> of  Specific  Activity  (From  Radial  to  Constant)  for 
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is  the  SOQu  curve  which  is  the  highest.  At  30  kt,  the  same  insensi- 
tivity to  the  choice  of  the  crossover  point  is  apparent  as  indicated 
in  Figures  6.8  and  6.0. 

6. 2 COMPARISON  OF  FREILING  MODEL  WITH  SIMPLE 

DISTRIBUTIONS 

There  is  a substantial  amount  of  uncertainty  associated  with 
the  distribution  of  fission  activity  with  particle  size.  As  discussed 
in  the  previous  section,  for  example,  the  Freiling  model  has  to  be 
modified  because  it  does  not  agree  with  the  observed  independence 
of  particle  size  above  a few  hundred  microns  insofar  as  specific 
activity  is  concerned.  In  addition  to  the  thermodynamic  effects 
involved,  there  must  also  be  a separate  mechanism  that  is  not 
modeled,  which  results  in  the  collection  of  fission  debris  by  the 
dust  with  a size -dependent  efficiency  factor. 

As  a further  investigation  of  the  sensitivity  of  the  predicted 
dose  to  model  assumptions,  modifications  were  made  to  the  PAM 
so  that  the  radial  component  of  the  activity  distribution  was  first 
purely  volatile  (surface  distributed)  and  second  purely  refractory. 

In  the  latter  case,  then,  the  entire  distribution  would  be  refractory. 
Figure  6. 10  shows  the  integrated  dose  along  the  hotline  at  3 Mt  for 
short  ranges  while  Figure  6. 11  is  for  intermediate  ranges.  The 
three  curves  shown  are  for  the  standard  DELFIC  calculation  and 
for  these  two  simpler  assumptions.  Predictions  of  the  30  kt  hot- 
line dose  are  shown  in  Figures  6. 12  and  6. 13.  At  both  yields,  the 
assumption  that  the  specific  activity  is  independent  of  particle  size 
(i.e. , all  of  the  decay  chains  are  purely  refractory)  is  on  the  aver- 
age 20  percent  higher  while  taking  the  radial  component  to  be 
volatile  (and  therefore  the  specific  activity  varies  with  particle 
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Figure  6.8.  Sensitivity  of  Hotline  Dose  to  Crossover  of  Specific 
Activity  (From  Radial  to  Constant)  for  30  kt  - 
Short  Range 
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Figure  6. 10.  Comparison  of  Freiling  Model  with  Assumptions 
of  Volatile  and  Refractory  Radial  Components  for 
3 Mt  Hotline  Dose  - Short  Range 
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Solid  Curve:  Modified  Freiling  Model 

Dashed  Curve:  Radial  Component  is 
Pure  Volatile 

Dot -Dashed  Curve:  Pure  Refractory 
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Figure  6. 11.  Comparison  of  Freiling  Model  with  Assumptions 
of  Volatile  and  Refractory  Radial  Components  for 
3 Mt  - Intermediate  Range 
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Figure  6. 13.  Comparison  of  Freiling  Model  with  Assumptions  of 
Volatile  and  Refractory  Components  for  30  kt 
Hotline  Dose  - Intermediate  Range 
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diameter  (d)  as  d~*  results  in  an  underestimate  of  the  dose  by  $0 
percent.  The  conclusion  to  be  drawn  here  is  that  for  short  and 
Intermediate  ranges,  the  simple  assumption  that  all  mass  chains 
are  purely  refractory  reproduces  satisfactorily  the  more  complex 
Freiling  model.  It  is  quite  possible,  however,  that  this  agreement 
is  significantly  influenced  by  the  choice  of  the  particle  number  dis- 
tribution, which  for  this  study  was  a power-law  (Section  2). 


6.3  REFERENCES 

1.  R.C.  Tompkins,  "Department  of  Defense  Land  Fallout 
Prediction  System,  Volume  V,  Particle  Activity,  ” 
DASA-1800-V,  U.S.  Army  Nuclear  Defense  Labora- 
tory, Edgewood  Arsenal,  Maryland,  1968. 
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Section  7 


SENSITIVITY  TO  UNCERTAINTIES 
IN  INDUCED  SOIL  ACTIVATION 

Even  though  the  induced  activity  is  not  a significant  compo- 
nent of  the  fallout  unless  the  fission  fraction  is  less  than  perhaps  50 
percent,  it  is  one  of  the  most  difficult  to  model.  This  section  is 
devoted  to  a discussion  of  the  uncertainties  associated  with  the  model 
of  induced  activity  that  is  presently  employed  in  the  Particle  Activity 
Module. 

According  to  Reference  1,  the  three  basic  assumptions  are: 

a.  all  neutrons  entering  the  soil  are  thermalized,  and 
then  captured 

b.  only  those  neutrons  that  are  seen  by  the  apparent 
crater  contribute 

c.  all  significant  soil  components  are  refractory  in  the 
fractionation  scheme. 

In  addition,  only  NTS  soil  is  provided  for. 

7. 1 COMPARISON  WITH  FISSION  ACTIVITY 

For  a surface  burst,  DELFIC  assumes  17  percent  of  the 
emitted  neutrons  are  captured  independent  of  yield  and  neutron  energy 
spectrum.  For  comparison,  we  have  the  work  of  Lessler  and  Guy 
(Reference  2).  Using  a Monte  Carlo  technique,  they  were  able  to 
determine  the  fraction  of  neutrons  captured  versus  radius  when  a 
mono-energetic  neutron  source  was  placed  on  the  surface  of  dry  NTS 
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soil.  A not -too -surprising  result  of  their  study  was  that  high-  and 
low-energy  neutrons  behave  differently.  Table  7. 1 lists  the  cumulative 
fraction  of  neutrons  that  are  captured  versus  distance  (along  the 
surface)  from  the  source  and  Figure  7. 1 is  a schematic  of  the  zones 
used  in  the  calculation.  Practically  all  of  the  neutrons  which  were 
captured  outside  the  5 m zone  were  transported  indirectly  through 
the  air. 


Table  7. 1 Fraction  of  Neutrons  Captured 


Distance 
From  Source 
(m) 

Cumulative  Fraction 

14  Mev 

0. 5 Mev 

5 

0.256 

0. 101 

15 

0.267 

0. 121 

30 

0.279 

0. 141 

45 

0.292 

0. 160 

900 

0.461 

0.386 

These  numbers,  which  include  the  effects  of  secondary  scattering  in 
air,  were  determined  assuming  ambient  air  and  soil  temperatures. 

If  there  is  pre-heating,  then  the  fractions  would  be  lower. 

(3) 

Table  7. 2 lists  the  radius  of  the  apparent  crater  for 
various  yields  and  the  fraction  of  neutrons  captured  within  the  crater 
boundary  for  surface  bursts.  These  fractions  were  obtained  by 
interpolation  in  Table  7. 1 which  was  derived  for  dry  NTS  soil  and 
the  assumption  is  made  in  this  analysis  that  the  values  derived  are 
valid  for  other  soils.  Note  the  energy  dependence  on  the  fraction 


t 


captured,  particularly  at  low  yields. 


Studies  of  the  soil  burden  in  the  stabilized  cloud  of  a sur- 
face burst  have  shown  it  to  be  only  a fraction  of  the  mass  removed  in 
the  cratering  process.  Consequently  the  fraction  of  captured  neutrons 
that  contribute  to  the  fallout  activity  is  going  to  be  less  than  indicated 

in  Table  7. 3.  For  example,  at  1 Mt,  the  soil  burden  of  the  stabilized 

(4) 

cloud  is  on  the  order  of  0. 3 Mt'  while  the  apparent  crater  mass  is 

about  5 Mt.  While  no  model  presently  exists  for  the  overall  cratering 

(5) 

and  entrainment  process,  it  has  been  observed'  that  the  soil  that  is 
ejected  and  not  entrained  tends  to  originate  away  from  the  burst  point. 
On  the  basis  of  this,  a simple  assumption  about  the  original  position 
of  the  entrained  soil  is  that  it  occupied  a volume  centered  on  the  burst 
point,  which  when  multiplied  by  an  average  bulk  density  for  soil  of 
2. 5 gm/cc  equals  the  stabilized  cloud  loading.  A further  assumption 
is  made  that  the  region  is  a paraboloid  with  depth  equal  to  half  the 
radius.  Table  7. 3 lists  the  fraction  of  neutrons  activating  the  soil 
in  the  stabilized  cloud.  It  is  these  fractions  that  are  to  be  used  in 
determining  the  contribution  of  neutron -induced  soil  activity  to  fall- 
out following  a surface  burst. 

Table  7. 3 Fraction  of  Neutrons  Activating  the  Soil  Burden 
of  the  Stabilized  Cloud 


w 

(Mt) 

Fraction 

14  Mev 

0.  5 Mev 

.01 

0.26 

0.11 

.10 

0.28 

0.13 

1.0 

0. 30 

0.16 

10. 

0.33 

0.21 
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The  use  of  these  fractions  versus  those  in  Table  7.2  results  in  only 
a 10  to  20  percent  decrease  in  the  estimate  of  activated  soil  in  the 
cloud.  However,  as  will  be  discussed  below,  the  original  location 
of  the  activated  soil  may  have  important  consequences  for  its  size 
distribution. 

As  an  example  of  the  relative  contributions  of  fission  debris 
and  activated  soil,  consider  a 300  kt  surface  burst  with  50  percent 

/g\ 

fission.  Using  hypothetical  neutron  output  spectra'  ' (in  which  the 

23 

fission  spectrum  is  normalized  to  a total  of  2.0  x 10  neutronsAt) 
the  number  of  neutrons  per  kt  emitted  above  and  below  4 Mev  is 

22  9* 

2. 64  x 10  and  2. 32  x 10  respectively.  This  partition  is  justified, 
since  below  4 Mev,  capture  cross  sections  are  roughly  energy  inde- 
pendent. In  this  analysis,  we  use  the  14  Mev  fractions  for  the 
neutrons  above  4 Mev  qnd  the  0.  5 Mev  fractions  for  the  remainder 
of  the  spectrum.  The  fracUon  in  each  energy  bin  which  activate 
entrained  soil  can  be  found  by  interpolation  in  Table  7. 3',  the  values 
are  29  percent  in  the  high  energy  bin  and  14  percent  in  the  low  energy 
bin.  Figure  7. 2 shows  the  dose  that  would  result  from  the  fission 

debris  and  the  high  and  low  energy  neutrons.  The  ordinate  is  a 

-1  2 

normalized  dose  in  units  of  R kt  mi  where  the  dose  rate  was  inte- 
grated from  time  of  arrival  to  H+48  hours.  To  facilitate  the  com- 
parison, it  was  assumed  first  that  the  size  distribution  of  the  neutron- 
induced  activity  was  the  same  as  that  for  the  fission  debris  and 
second  that  the  fission  decay  rate  was  proportional  to  T ' where 
T is  time  in  hours  after  burst.  This  calculation  was  done  for  dry 
NTS  soil.  Reference  1 provides  the  decay  rate  per  mole  of  neutrons 
(14  Mev  and  0.  5 Mev)  both  individually  for  each  soil  element  and 
collectively.  Thus  the  curves  in  Figure  7. 2 required  only  the  number 
of  neutrons  in  each  energy  bin.  Note  that  although  the  number  of  high 
energy  neutrons  emitted  is  only  10  percent  of  the  total,  they  contribute 
over  20  percent  of  the  total  induced  activity. 
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Figure  7. 2 Comparison  of  Fission  and 

Neutron-Induced  Soil  Activity 
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The  foregoing  analysis  has  not  been  based  on  neutron  energy 
spectra  from  actual  weapons,  hence  caution  should  be  exercised  in 
attempting  to  draw  generalized  conclusions. 

7. 2 SOIL  VARIATION  EFFECT  ON  INDUCED  ACTIVITY 

The  Particle  Activity  Module  is  presently  set  up  only  to 
handle  NTS  soil.  However,  the  chemical  composition  is  substantially 
different  than  that  reported  for  NTS  soil  in  Lessler  and  Guy’s  report 
(Reference  2).  The  DELFIC  profile  is  much  higher  in  hydrogen 
(with  an  atomic  fraction  of  0.42  versus  0. 16  in  Reference  2)  and  con- 
versely a lower  proportion  of  Na  and  Al.  The  discussion  to  follow 
assumes  the  NTS  profile  of  Reference  2. 

In  addition  to  the  NTS  soil,  four  average  profiles^  have 
been  considered  and  are  the  following: 

a.  Chernozem  — an  average  of  five  samples  of  this 
soil  group  from  the  USSR 

b.  Podzol  — an  average  of  two  samples  of  this  soil 
group  from  the  USSR 

c.  USSR  — average  of  all  samples  in  USSR 

d.  PRC  — average  of  five  Peoples  Republic  of  China 
samples. 

The  data  given  in  Reference  5 are  only  for  certain  elements  and  are 
in  terms  of  weight  fractions  and  are  listed  in  Table  7.4.  In  particular 
the  water  fraction  is  omitted.  In  order  to  make  a comparison  with 
the  NTS  calculation  of  Reference  2,  the  following  has  been  assumed: 

a.  the  same  weight  fraction  (0. 06)  of  water  as  for  dry 
NTS  soil 
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b.  the  same  weight  fraction  (0.003)  for  other  hydrogen 

c.  a mean  atomic  weight  determined  from  the  data  of 
Reference  2 for  NTS,  Chester  and  Dade  soils. 


Table  7. 4.  Weight  Fractions  for  Representative  Soils 


Element 

Chernozem 

Podzol 

USSR 

PRC 

NTS 

Na 

■ ■ 

0.0188 

0.0116 

pph 

0.0134 

A1 

■ ■ 

0.0842 

Si 

1 

0.3084 

! 0.3036 

0.27561 

K 

’ TT»l§t. 

0.0200 

M'WtYr  nt 

mm. 

Ti 

0 

0 

0.0058 

Mn 

iKoi 1 * 

0.0011 

0.0004 

Fe 

0.0313 

Mg 

KVjfjfVI 

Ca 

0.0293  | 

0.0144  | 

0.0220  1 

warn 

0.0866 

This  latter  factor  is  determined  as  follows: 


i 

I — , __FW. 

| A = 0.94  2^-j-i-  +0.06/18 

i i ** 


where  Aj  is  the  atomic  weight  of  the  ith  element  and  FW.  is  its  weight 

fraction.  The  second  term  represents  the  water  contribution  and  the 

summation  in  the  first  term  is  renormalized  since  the  weight  fractions 

were  determined  with  the  water  removed.  The  mean  value  of  A”1 

-2 

was  found  to  be  5. 21  x 10  . Using  this,  the  atomic  fractions  were 
determined  as  follows: 
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FAt  = 0.94  FW^A-1  x 

and  are  listed  in  Table  7.  5.  Figure  7. 3 compares  the  ratio  of  the 
induced  activity  for  each  soil  relative  to  the  fission  activity  for 
a 300  kt  surface  burst  with  50  percent  fission  fraction.  The  fission 
and  fusion  neutron  spectra  are  those  of  EM-1  and  are  not  intended  to 
be  typical  of  modern  weapons.  The  activity  is  an  integral  of  the 
dose  rate  from  time  of  arrival  to  H+48  hrs.  In  determining  the  rela- 
tive contribution  made  by  the  activated  soil,  it  was  assumed  to  be 
constant  over  the  dust  particle  size  spectrum.  The  Podzol  soil  has 
the  largest  effect  while  the  PRC  soil  is  lowest  — being  only  20  percent 
of  the  Podzol-induced  activity.  Thus,  there  is  a factor  of  five  varia- 
tion but  even  the  largest  is  less  than  20  percent  of  the  fission  activity 
for  this  particular  case  of  50  percent  fission  fraction. 


Table  7. 5.  Atomic  Fractions  for  Representative  Soils 


Soil  Type' 


Chernozem 

Podzol 

USSR 

PRC 

Dry  NTS 

Na 

0.00675 

0.00910 

0.00230 

0.0105 

A1 

0.0466 

0.0514 

0.0577 

0.0510 

Si 

0. 193 

ran 

0.196 

0.201 

0.177 

K 

0.00810 

0. 00925 

0.00745 

0.00320 

0.0137 

Ti 

0 

0 

0.000075 

0.00218 

0.0009 

Mn 

0.000066 

0.00036 

0.00020 

0.00013 

0. 000098 

Fe 

0.0108 

0.0101 

0.0113 

0.0157 

0.0036 

Mg 

0.00677 

0.00068 

0.0055 

0.00346 

0.0033 

Ca 

0.0132 

0.0065 

0.0099 

0.0013 

0.0390 

An  important  factor  in  assessing  the  importance  of  induced 
activity  to  fallout  is  the  percent  of  water  in  the  soil.  As  pointed  out 
in  Reference  2,  hydrogen  has  a capture  cross  section  that  is  approxi- 
mately three  times  higher  than  that  for  the  other  soil  constituents. 
Consequently,  it  serves  to  diminish  the  number  of  neutrons  that  would 
activate  other  elements.  In  the  preceding  comparison  among  various 
soils,  it  was  assumed  that  each  had  the  same  fraction  of  water  (6  per- 
cent) as  in  the  dry  NTS  soil  used  in  Reference  2.  Thus,  in  estimating 
the  number  of  activated  atoms,  all  that  was  necessary  was  to  multiply 
the  values  supplied  in  Reference  2 using  dry  NTS  by  the  ratio  of 
atomic  fractions  of  a given  soil  to  that  for  the  dry  NTS  soil. 

7.3  VARIABLE  WATER  CONTENT 

A comparison  of  wet  versus  dry  NTS  soil  by  Lessler  and 
Guy  Indicated  that  an  increase  of  the  water  content  from  6 to  17 
percent  by  weight  decreased  the  activation  by  35  percent.  The 
following  shows  the  mathematics  of  this  effect  and  indicates  how  a 
variable  water  content  would  affect  the  soil  activation.  The  base- 
line is  water  content  of  6 percent  by  weight. 

It  is  first  necessary  to  determine  the  dependence  of  the 
hydrogen  atomic  fraction  due  to  water.  Following  the  above  methodol- 
ogy, the  mean  atomic  weight  inverse  is  given  by: 


— r *-w. 

A * (l  - FWH.O>  2 -IT  * ™H. 


,0/18 


where  FW„  is  the  weight  fraction  of  water  and  the  summation  is 

h2u  — n 

over  the  elements  listed  in  Table  7.4.  For  the  baseline  case,  A 

-2 

was  found  to  be  5. 21  x 10  which  implies  that  the  summation  in  the 

-2 

above  expression  is  5. 19  x 10  . Thus,  for  other  water  weight  fractions, 
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A"1  = (1  - FWh  q)  (5. 19  x 10“*)  + FWH  0/18. 

2 2 


The  atomic  fraction  of  water  (HgO)  is  then 


fah20  - fwh2o/(a  * “>• 

Consequently  the  atomic  fraction  of  the  hydrogen  in  the  water  is 


FAjj  = 2/3(FAh^0)  = (2/3)FWh^0/(A‘x  18) 

where  A”^  is  dependent  on  FW^q.  For  the  baseline  case,  FA^  is 
.043,  while  for  wet  NTS  soil  it  is  equal  to  . 120.  If  the  water  con- 
tent by  weight  were  30  percent,  FA^  would  be  . 209. 

The  total  number  of  activated  soil  atoms  can  be  expressed 
as: 


fc~l-*-rAH 


where  x is  the  ratio  of  the  absorption  cross-section  by  hydrogen  to 
the  activation  cross-sections  of  the  other  soil  elements.  To  deter- 
mine x,  the  fact  that  wet  versus  dry  NTS  soil  reduces  f by  35  per- 
cent is  used,  i.  e. , 


fWET  = 0.65fDRV 
c c 


or 


0.65 


1 - . 120  x 

1 - . 643  x 


which  results  in  x being  equal  to  3.  8.  For  any  other  water  content, 
then,  the  decrease  in  activation  from  the  baseline  case  is  given  by: 
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1-3.8  FA, 


•JDffT 


= 1.20(1  -3.8  FAjj) 


For  water  contents  of  10,  20  and  30  percent  by  weight,  this  fraction 
is  0. 88,  0. 56  and  0. 25  respectively.  Figure  7. 4 shows  the  dependence 
on  water  content  using  Podzol  soil. 

7. 4 HEIGHT-OF-BU RST  EFFECTS  ON  INDUCED  ACTIVITY 

The  Particle  Activity  Module  accounts  for  HOB  variations 
by  calculating  the  fraction  of  total  solid  angle  intersected  by  the 
apparent  crater.  Relative  to  a surface  burst,  the  dependence  is 
given  by 


F = l - 8 ^ 

V*.24  HOB  ^ - 234  HOB  + 4225 
s s 

0 < HOB  < 36  ft 
s 

= 0,  HOB  > 36  ft 
s 

where  HOB  is  the  actual  height  of  burst  scaled  to  1 kt  by  the  con- 
ventional  cube-root  of  yield  scaling.  However,  an  praminatinn  of 
the  HOB  dependence  in  Reference  3 of  apparent  crater  dimensions 
indicated  the  log  of  the  apparent  crater  volume  to  fall  off  linearly 
Wth  HOBg  up  to  10  feet  with  the  volume  at  HOBg  = 10  ft  being  0.01 
of  the  crater  volume  from  a surface  burst.  If  it  is  assumed  that 
both  the  crater  radius  and  its  depth  have  the  same  exponential 
dependence  on  HOB,  then 


143 


Ra(HOB  ) = Ra(HOB  = 0)  • exp  (-HOB  /8.«) 

8 S 8 


where  Ra  is  apparent  crater  radius  and  all  dimensions  are  in  feet. 
Furthermore,  above  a scaled  height  of  burst  of  approximately  15  feet 
the  crater  formation  mechanism  changes  from  ejection  to  compac- 
tion. Thus  the  HOB  dependence  of  the  induced  activity  should  go  to 

zero  at  HOB  = 15  feet.  Using  the  same  rationale  as  presently  con- 
8 

tained  in  the  PAM,  the  formulation  of  F should  be 


F - 1. 


HOB 

B 

VhOBJ*  + (Ra  W1/3)2 

D 


0 < HOB  < 15  feet 

8 

F = 0,  HOB  > 15  feet 
s 

where  Ra  has  the  dependence  on  HOB  given  above  and  W is  the 

s 

yield  in  kt.  Figure  7.  5 compares  the  DELFIC  formulation  with 
this  calculation  for  a 1 Mt  yield.  It  is  apparent  that  DELFIC  is 
presently  overestimating  the  correction  factor  above  a few  scaled 
feet. 

Not  considered  in  the  above  analysis  is  the  induced  activity 

in  the  surface  layer  of  dust  part  of  which  will  be  entrained  even  at 

HOB„  for  which  no  crater  is  formed.  However,  this  dust  was 
a 

originally  located  down  to  a depth  of  only  a few  millimeters.  A sim- 
ple calculation  of  the  transmission  of  neutrons  through  soil  indicates 
that  less  than  five  percent  are  captured  in  the  first  5 mm.  Conse- 
quently, the  induced  activity  in  the  surface  layer  has  not  been  con- 
sidered. However  further  calculations  should  be  performed. 


i 
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8 16  24  32  36  40 

Scaled  HOB  (Ft) 


Figure  7.  5 HOB  Correction  Factor  for  Induced  Activity 
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7.  5 SIZE  DISTRIBUTION  EFFECTS  ON  INDUCED  ACTIVITY 

DELFIC  assumes  that  the  activated  soil  is  refractory, 
meaning  that  it  is  volume  distributed  throughout  the  fallout  material. 
However,  as  has  been  pointed  out  in  Section  7. 1,  most  of  the  neutrons 
are  captured  close  to  the  burst  point.  Crater  material  that  experiences 
a peak  stress  greater  than  a megabar  can  be  expected  to  be  vaporized, 
while  outside  this  region  the  dominant  effect  is  to  break  up  the  material 
mechanically.  For  a 1 Mt  surface  burst  on  rock,  the  1 megabar 
point  occurs  at  a range  of  about  13  meters  within  which  occur  80  to 
90  percent  of  the  neutron  captures  that  contribute  to  fallout.  The 
fraction  of  entrained  material  originating  in  this  region  is  only  on 
the  order  of  2 percent.  The  remainder  may  be  subjected  to  some 
degree  of  melting  as  it  is  ejected  into  the  fireball  but  there  is  probably 
little  vaporization  occurring  in  the  fireball.  Consequently,  it  would 
appear  that  a reasonable  alternative  to  the  DELFIC  assumption  is 
that  the  activated  soil  is  volatile,  i.  e. , it  is  preferentially  area 
distributed  on  the  fallout  particles.  To  determine  the  sensitivity  to 
this  assumption,  the  Particle  Activity  Module  was  modified  in  the 
following  manner  for  the  3 Mt  case: 

a.  the  fission  activity  was  set  to  zero 

b.  the  induced  activity  was  changed  to  a volume  distri- 
bution with  size  below  lOOp  and  area  distribution 
above  100p  . 

Figure  7.6  shows  the  comparison  with  the  current  DELFIC  assumption. 
The  curves  are  for  integrated  dose  versus  range  along  the  hotline  as 
defined  in  Section  2.  At  close-in  distances,  the  DELFIC  prediction 
for  the  contribution  of  induced  activity  is  higher  as  would  be  expected 
since  *he  assumption  that  the  activated  soil  can  be  treated  as 


4# 


Figure  7.6  Comparison  of  Alternative  Fractionation 
Assumptions  for  Induced  Activity  at 
Close-In  Ranges 


148 


r,;.  ■ 


refractory  places  a higher  fraction  of  this  component  of  the  total 
activity  on  the  larger  particles.  At  zero  range,  the  DELFIC  pre- 
diction is  an  order  of  magnitude  higher  with  the  two  crossing  at 
about  25  km.  The  total  activity  is  also  shown  on  Figure  7.3.  At 
ranges  beyond  25  km,  the  model  described  above  results  in  a higher 
prediction  as  shown  with  less  resolution  on  range  in  Figure  7. 7. 
However,  at  ranges  greater  than  100  km,  the  induced  activity  contri 
button  predicted  either  way  is  below  10  R for  this  particular  weapon 
configuration. 
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Appendix  A 

USE  OF  THE  STABILIZED  CLOUD  TAPE  FROM 
SUBROUTINE  RSXP  IN  THE  CLOUD  RISE  MODULE 


i 

i 

i 

4 

« 

1 


According  to  DASA-1800-IH  (Revised),  the  RSXP  debug 
printout  gives  properties  of  the  cloud  wafers  before  they  are 
sectioned  in  the  horizontal  plane.  The  printout  column  headings  and 
their  definitions  are  given  in  Table  A.  1. 


Table  A.  1.  Printout  Explanation 


TIM 

ALT 

RAD 

DIAM 

MASS 

DZ 

ZLOW 

VOL 

MBT 


Time  (sec) 

Altitude  of  wafer  center  of  mass  (m  above  msl) 
Radius  (m) 

Particle  size  class  midrange  diameter  &/m) 

Total  particulate  mass  in  the  wafer  (kg) 

Wafer  thickness  (m) 

Wafer  bottom  altitude  (m  above  msl) 

Wafer  volume  (m  ) 

(Always  = 1)  signifies  that  both  wafer  top  and  bottom 
have  been  processed 


IF  LAG 


A parameter  that  signifies  whether  a wafer  is  part 
of  the  cloud  cap  or  stem.  If  it  is  totally  or  partially 
in  the  stem,  further  vertically  subdivided  wafers 
are  printed  out  next. 

IFLAG  = 1 no  further  subdivision  required 


IFLAG  = 2 further  subdivision  required 


A typical  example  of  an  RSXP  printout  is  the  section  of  the 
3 Mt  data  shown  as  Table  A.  2.  First,  note  that  all  the  volumes  are 
0.000.  Secondly,  the  height  of  the  cloud  at  stabilization  is  given  as 
2.2098+04  m,  which  means  that  ALT  is  the  wafer  top  height,  not  the 
wafer  middle  height  since  DZ  is  given  as  2.059+02  and  ALT  is  given 
as  2.210+04  in  the  last  row  of  the  table.  Thirdly,  and  most  impor- 
tantly, the  values  of  DZ,  ZLOW,  and  VOL  seem  to  be  retarded  one 
row  in  the  table.  It  appears  the  dotted  line  shown  follows  the  correct 
values.  Using  the  values  outlined, 

ALT  - ZLOW  = 2. 190+04  - 2. 169+04  = 0.021+04  = DZ. 

This  is  more  apparent  in  prior  sections  of  the  printout  where  DZ 
varies  over  the  particle  size  class.  And  finally,  the  particle  size 
listed  is  the  lower  boundary  diameter  of  each  size  class,  not  the 
midrange  diameter. 


Table  A.  2.  RSXP  Debug  Printout  for  3 Mt 


f\.<\'V?VfvfUf\.rv<vAjr\.cj<vn.»vr.<,».<\.c  f\i  «v  A *■  «\.  f\.  fur* *.  rv, 
M 4BOOOOOOOOOOOOOOOOOOCOOOaOOOOC>OOOOOOOOC  O o c c c c o OO  o 

^ ecc,“ccc,ecce^ccc-c*ccc'e-“rc-<rccirccrc‘r<rcc-ccccccccec  r V*c  c 
. *cxrrxcixx=.a  a ~ z.  z i s - : ~ 

H rr?3j^3j  j'r'a  jr 

crrcceccrrcc^cecrcrrrcervcrcccrcccccrcrcccrcc  «l«l«  • 


entewi'io  wnosrr 


Appendix  B 

DISCREPANCIES  IN  THE  GROUNDED  PARTICLES  TAPE 

According  to  DASA-1800-IV,  the  grounded  particles  tape 
contains  the  X and  Y coordinates  of  each  grounded  particle,  time  of 
impact,  particle  diameter,  mass  per  unit  area,  and  an  area  factor 
associated  with  the  particle.  The  wind  hodograph  used  with  each 
yield  is  shown  as  Table  B.  1. 

Table  B.l.  Wind  Field  Data 


I 

M) 

t-  ■ ' 

f 

? A planar  topography  was  assumed,  and  there  were  no  up- 

ft 

? drafts  or  other  localized  wind  phenomena  input.  Consequently,  an 

increasing  value  of  X and  a decreasing  value  of  Y correspond  to  an 
increasing  time  of  impact.  Tables  B.  2 and  B.  3 are  reproductions 


Vector  Altitude,  ZV(J) 

VX(J) 

VY(J) 

0.00000 

4. 55700f00 

-8. 19000-01 

3.04800+03 

4.557004-00 

-8. 19000-01 

6. 09600+03 

6. 758004-00 

-2.458004-00 

9. 14400+03 

1.008604-01 

-1. 792004-00 

1. 21920+04 

1.203204-01 

-4.403004-00 

1. 52400*04 

1. 203204-01 

-4. 40300+00 

1. 82880f 04 

1.510404-01 

-2. 66200+00 

2. 133604-04 

1.740804-01 

-3. 12300+00 

2. 438404-04 

2.048004-01 

0.00000 

2.743204-04 

2.304004-01 

0.00000 

3.048004-04 

2.560004-01 

0.00000 

of  a listing  of  the  3 Mt  grounded  particles  tape.  Various  rows  of 
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leB.2.  Section  of  3 Mt  Tape  IPOUT 
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ccrrccccccccccccccccccccpcsceccccccccccrrcrccccccccccrcpc 

ecec  c oc  ccccccc  c c c cccc  c c c c c c ec  ccccececc  c cee  — — — — ——  ——  — — — — — — — 

ex  ceccccccccccccccccccccccccccccccccccc  c c c p p»  p p p p p p p p hP  ^ h f 

+ r cv  <\  <v  f c.  «\  r.  n p <\  r.  r r r 
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Table  B.  3.  Section  of  3 Mt  Tape  IPOUT 
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data  are  marked  with  symbols  (*,+  , +)  to  Illustrate  pairs  of  data 
which  are  inconsistent.  Each  pair  is  for  a different  particle  size, 
and  yet  each  _ it  has  the  same  error.  For  an  increasing  X and 
decreasing  Y,  T is  decreasing  instead  of  increasing. 

Also,  the  arrow  symbol  ( ) on  Table  B.  3 shows  where  the 
X,  Y,  and  T coordinates  all  increase  by  a factor  of  10  for  no  obviously 
apparent  reason. 

And  finally,  also  on  Table  B.  3,  the  symbol  o shows  two 
rows  of  data  which  not  only  exhibit  the  increasing  X-decreasing 
Y -decreasing  T error,  but  also  shows  a wafer  with  a larger  radius 
hitting  the  ground  before  a wafer  with  a smaller  radius.  By  the 

v 

logic  of  the  Cloud  Rise  Module,  this  is  impossible. 

It  appears  that  there  is  an  error  associated  with  the  time  of 
impact  calculations.  Similar  errors  were  discovered  in  the  10  Mt 
data.  These  errors  were  not  observed  in  the  30  kt  data,  nor  were 
they  observed  in  the  3 Mt  or  10  Mt  data  for  wafer  radii  multipliers 
much  less  than  886.  The  30  kt  data  contain  very  few  wafer  radii 
multipliers  of  886.  Either  the  time  error  is  yield  dependent,  or 
wafer  radii  dependent. 
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Appendix  C 

FORTRAN  LISTING  OF  PROGRAM  L INTAP 
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